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The  radial  structure  of  a diffusion  controlled  oxygen  discharge 
at  equilibrium  is  mathematically  modeled,  and  is  solved  analytically 
by  assuming  a Maxwellian  electron  energy  distribution  and  then  using 
a Bessel  substitution.  A model  consisting  of  electrons.  O',  and  0^ 
ions  in  background  O2  at  300®K  along  with  the  attendant  reactions, 
dissociative  attachment  (O2  + e"  ->■  0“  +0),  its  reverse  (O'  + 0 -»• 

O2  + e'),  and  ionization  (O2  + e'  0^  + 2e'),  is  solved  analytically 
as  a basic  characteristic  value  problem  for  the  electron  temperature. 
It  is  found  that  for  a tube  radius  of  1.1  cm  and  at  pressures  of  0.2, 
0.5,  1.0,  and  5.0  torr,  the  equilibrium  electron  temperatures  are 
3.7  eV,  2.2  eV,  1.6  eV,  and  1.0  eV,  respectively.  The  results  also 
indicate  that  the  O'  ion  is  the  numerically  dominant  negative  charge 
carrier  and  that  it  diffuses  radially  inward.  The  analytic  technique 
is  further  applied  to  a four  species,  four  reaction  discharge  model 
which  contains  the  0^  ion  and  the  charge  transfer  reaction 
(O2  + O'  -►  O2  + 0).  An  analytic  solution  is  found  for  this  model. 

The  pressure-radius  product,  pR,  is  not  a scaling  parameter.  The 
solution  of  these  models  serves  to  prove  the  viability  of  a simple 
analytic  scheme  which  may  be  used  to  solve  other  discharge  models 
which  contain  more  than  four  charged  species  and  more  than  four 
accompanying  reactions. 
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RADIAL  STRUCTURE  IN  THE  LOW  PRESSURE 
EQUILIBRIUM  OXYGEN  DISCHARGE 

I . Introduction 

The  universe  is  almost  entirely  composed  of  plasmas.  Stars  are 
gigantic  plasmas.  The  relative  void  between  the  stars  is  filled  with 
low  density  plasmas.  Yet,  man  has  evinced  little  interest  in  this 
ubiquitous  phenomenon  until  a comparatively  short  time  ago.  Today, 
plasma  research  has  blossomed  into  one  of  the  dominant  fields  of 
scientific  inquiry. 

A plasma  is  an  electrically  conducting  medium  composed  in  general 
of  charged  particles,  both  positive  and  negative,  whose  properties 
depend  on  the  collective  behavior  of  the  particles,  as  distinct  from 
the  individual.  The  present  day  theoretical  effort  to  advance  man's 
knowledge  of  the  plasma  state  is  founded  on  Lorentz's  work  in  applying 
Boltzmann's  particle  kinetic  theory  to  the  behavior  of  electrons  in 
metals  and  on  wave  particle  kinetics.  Since  the  1950's,  space  ex- 
ploration, electronic  devices,  gas  discharge  lasers,  and  the  desire 
to  achieve  a controlled  thermonuclear  reactor  have  all  stimulated 
interest  in  plasmas. 

Some  common  examples  of  plasmas  are  auroras,  lightning,  welding 
arcs,  and  fluorescent  lighting.  The  last  is  an  example  of  an  ex- 
tended positive  column  discharge  in  an  insulated  glass  tube.  The 
discharge  tube  is  an  easy  and  convenient  means  of  producing  a plasma. 
An  electrical  potential  is  applied  across  a tube  filled  with  gas  and 
the  electrons  are  driven  to  high  kinetic  energies.  These  energetic 
electrons  in  turn  may  cause  ionization  of  the  gas  particles;  thus  a 
plasma  is  born. 
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Eventually,  an  equilibrium  may  be  establiihed  within  the  discharge 
such  that  the  production  and  loss  of  charged  particles  achieves  a 
balance.  Charged  particles  are  lost  and  reborn  through  the  reactions 
in  which  they  participate.  However,  particles  may  also  be  lost  as  a 
result  of  radial  migration  to  the  walls  of  the  discharge  tube.  In  a 
plasma  the  different  mobilities  of  electrons,  positive  ions,  and 
negative  ions  induce  electric  fields  which  tend  to  retard  the  migration 
of  the  electrons  and  negative  ions  and  which  tend  to  accelerate  the 
migration  of  the  positive  ions  toward  the  tube  wall.  In  this  special 
case  of  radial  migration,  the  combined  diffusion  of  the  electrons  and 
ions  forced  by  the  space  charge  electric  field  is  called  "ambipolar 
diffusion"  (Ref  18:417). 

An  example  of  a discharge  system  which  has  been  observed  experi- 
mentally off  and  on  over  the  past  fifty  years  and  which  is  currently 
being  investigated,  is  the  oxygen  discharge. 

Oxygen  Discharge  Models 

The  oxygen  discharge  contains  positive  ions,  negative  ions,  and 
electrons.  It  has  been  the  subject  of  considerable  research  over  the 
past  fifty  years  because  it  is  an  atmospheric  gas.  Recently,  it  has 
been  discovered  that  oxygen  can  be  made  to  lase.  This  development  has 
stimulated  experimentation  on  the  oxygen  discharge  system  in  an  attempt 
to  unravel  the  underlying  plasma  mechanisms,  and  thus  to  lay  the 
groundwork  for  the  design  of  an  efficient  lasing  system. 

The  oxygen  discharge  system  is  made  up  of  more  than  six  different 
species  and  more  than  170  processes  (Ref  2:263-267).  Because  the 
system  is  relatively  complex,  it  has  withstood  attempts  at  rigorous 
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mathematical  analysis.  However,  some  computer  models  have  been 
recently  developed  to  investigate  it  in  detail  (Ref  2). 

The  major  species  involved  in  the  oxygen  discharge  may  be  quali- 
tatively surmised  by  considering  the  electronegative  character  of 
oxygen  and  the  production  and  loss  mechanisms  of  the  various  species. 

At  a gas  temperature  of  around  SOCK  atomic  oxygen,  0,  has  an  electron 
affinity  of  about  1.5  eV;  and  molecular  oxygen,  O2,  has  an  electron 
affinity  of  approximately  0.5  eV.  Thus,  it  may  be  expected  that  both 
O'  and  O2  ions  will  be  formed  in  the  presence  of  electrons. 

Beginning  with  02*  electrons,  and  an  electric  field,  it  is  possible 
to  produce  O'  ions  via  dissociative  attachment  (0^  + e'  -+  O'  + 0) . The 
electron  energy  necessary  for  this  reaction  is  the  energy  required  to 
dissociate  molecular  oxygen  minus  the  electron  affinity  of  atomic 
oxygen  or  about  3.8  eV.  Once  O'  ions  are  produced,  a continuous  charge 
exchange  process  may  begin  (O'  + O2  ^ O2  + 0)  which  produces  O2  ions. 

Positive  ions  are  produced  by  ionization  of  both  O2  and  0.  Elec- 
trons with  energies  greater  than  about  12  eV  will  produce  0^  ions 
(O2  + e'  ->■  O2  + 2e'),  and  electrons  with  energies  greater  than  about 
13.5  eV  will  produce  0'*’  ions  (0  + e'  -►  0^  + 2e') . These  ions  may  also 
be  expected  to  participate  in  charge  transfer  reactions  (0^  + 0 ? 

O’*"  + 02).  Dissociative  ionization  reactions  require  electrons  with 
energies  greater  than  17  eV. 

Loss  of  oxygen  ions  can  occur  by  recombination.  Radiative  re- 
combination events  have  very  low  cross  sections  (Ref  20:521).  There- 
fore, the  only  important  loss  mechanisms  are  dissociative  recombination 
and  mutual  neutralizations  with  dissociation.  However,  Thompson  has 
shown  that  these  recombination  loss  rates  are  much  less  than  the 
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measured  wall  loss  for  diffusion  dominated  discharges  in  which  the 
volume  to  surface  ratio  is  much  less  than  one  (Ref  20:521). 

Loss  of  oxygen  ions  can  also  occur  by  detachment.  However,  many 
of  these  reactions  require  ion  temperatures  greater  than  1000°K; 
therefore,  they  are  not  significant  loss  mechanisms  for  discharges 
maintained  at  room  temperature. 

Of  the  species  involved  in  the  oxygen  discharge,  Thompson  be- 
lieves that  the  dominant  ions  are  0”  and  0^  negative  ions,  and  0^ 
positive  ions  (Ref  20:519).  Using  electrons,  and  the  0"  and  0^  ions, 
Sabadil  has  modeled  the  oxygen  discharge  utilizing  a numerical  scheme 
(Ref  17).  Recently,  Lee  has  investigated  an  analytic  solution  to 
discharge  systems  (Ref  10).  This  thesis  will  attempt  to  verify  that 
his  mathematical  technique  and  a limited  account  of  the  plasma 
chemistry  gives  reasonable  results  for  the  oxygen  discharge. 

The  first  simplified  model  analyzed  will  be  a three  species 
system  containing  electrons,  and  the  O"  and  0^  ions.  The  reactions 
considered  are  dissociative  attachment  (O2  + e”  -►  0”  + 0),  its  re- 
verse (O'  + 0 O2  + e'),  and  ionization  (O2  + e'  -►  0^  + 2e'). 

Second,  to  this  system  will  be  added  the  charge  transfer  reaction 
(O2  + O'  O2  + 0) , which  is  the  primary  means  for  producing  the  fourth 
species,  0^  (Ref  20:522).  The  third  and  fourth  models  will  include 
other  reactions  added  to  the  four  species  discharge. 

The  discharge  will  be  assumed  to  be  contained  within  a long 
cylindrical  tube  whose  dimensions  are  greater  than  the  Debye  length. 

The  Debye  length,  X^,  is  the  maximum  distance  in  a plasma  over  which 
a large  difference  in  the  number  of  positive  and  negative  charges  can 
occur. 
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Te  is  the  electron  temperature  in  °K,  and  Ne  is  the  number  density  of 

_3 

electrons  in  m . The  Debye  length,  is  given  in  meters.  For 
electron  temperatures  of  4 x 10^  “K  and  electron  densities  of  10^®  m"^, 
the  resulting  Debye  length  is  less  than  0.05  cm.  The  radius  of  the 
discharge  tube  in  this  thesis  is  1.1  cm  and  the  length  is  considered 
to  be  greater  than  1.0  m.  Therefore,  charge  quasi-neutrality  will  be 
assumed  to  exist  in  all  the  discharge  models  examined. 

Also,  it  will  be  assumed  that  electrons  in  the  discharge  may  be 
described  by  a Maxwellian  energy  distribution.  Actually,  Dettmer  pro- 
vides evidence  that  the  electron  energy  distribution  is  between  a 
Maxwellian  and  a Druyvesteyn  distribution  (Ref  2:221).  A Maxwellian 
distribution,  successful  in  past  discharge  models,  will  be  assumed 
here  (Ref  11). 

An  analytic  solution  will  be  obtained  for  both  the  three  and  four 
species  models. 

Sequence  of  Development 

In  Chapter  II,  discharge  equations  developed  from  the  first  three 
moments  of  the  Boltzmann  equation  for  a three  species  model  will  be 
solved  analytically  as  a basic  characteristic  value  problem  for  electron 
temperature.  It  will  be  assumed  that  the  electrons  are  Maxwellian, 
that  the  discharge  tube  has  cylindrical  syimetry,  that  there  is  no 
applied  magnetic  field,  and  that  self-electric  fields  are  negligible. 

In  order  to  obtain  numerical  results  from  the  analytic  solution 
of  the  three  species  model,  it  is  necessary  to  know  how  the  dissociative 
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attachment  rate  and  ionization  rate  depend  on  the  electron  temperature. 
Chapter  III  and  Chapter  IV  report  the  modeling  of  these  reaction  rates 
as  functions  of  the  electron  temperature. 

Chapter  V presents  the  numerical  results  from  the  solution  of  the 
three  species  oxygen  discharge  model  using  the  constants  pertinent  to 
oxygen  and  the  results  for  the  reaction  rates  obtained  in  Chapter  III 
and  Chapter  IV. 

In  Chapter  VI  the  oxygen  discharge  is  modeled  as  a four  species, 
four  reaction  system  and  solved  analytically,  again  as  a basic  charac- 
teristic value  problem.  Numerical  results  are  also  presented. 

In  Chapter  VII  some  other  models  are  brought  forth,  and  Chapter 
VIII  and  Chapter  IX  conclude  the  thesis  with  a discussion  of  the 
analytic  technique  used  to  evaluate  the  various  models  and  with  recom- 
mendations about  other  areas  that  could  be  appropriately  pursued. 

Notation 

All  equations  are  identified  by  sequence  number  according  to  the 
chapter  or  appendix  in  which  they  are  located,  e.g.  (2.11)  indicates 
equation  eleven  in  Chapter  II.  Figures  are  identified  in  the  same 
manner.  Tables  are  identified  simply  by  a Roman  numeral  in  the  order 
that  they  appear  in  the  thesis. 


II . An  Equilibrium  Model  for  the  Three  Species, 

Three  Reaction  Oxygen  Discharges 

A three  species,  three  reaction  model  of  the  O2  discharge  at 
equilibrium  can  be  developed  from  the  first  three  moments  of  the 
Boltzmann  equation.  By  specifying  a cylindrical  discharge  geometry 
and  by  utilizing  the  assumptions  relevant  to  the  given  discharge,  a 
system  of  differential  equations  which  models  the  discharge  may  be 
obtained.  These  differential  equations  will  be  called  "general  dis- 
charge equations." 

The  three  species  considered  are  electrons,  negative  atomic  oxygen 
ions,  and  positive  molecular  oxygen  ions.  They  are  represented  as 
e”,  0”,  and  Also  involved  are  the  background  neutral  oxygen 
molecules  O2,  and  the  neutral  oxygen  atoms  0. 

The  three  reactions  used  in  the  model  are  as  follows: 


1 


Name 

Reaction 

Rate (sec) 

Dissociative  Attachment 

CSJ 

0 

+ e"  0"  + 0 

ct(Tg) 

(2.1) 

Atomic  Detachment 

0“ 

+ 0 -*■  O2  + e" 

Y 

(2.2) 

Molecular  Ionization 

02 

+ e"  -►  O2  + 2e" 

3(Tg) 

(2.3) 

Both  the  dissociative  attachment  rate  a and  the  ionization  rate  B are 
modeled  as  functions  of  the  Maxwellian  electron  temperature  T^.  Also, 
a and  B are  linear  functions  of  the  background  gas  number  density  NO2. 
The  detachment  rate  y will  be  assumed  a function  of  the  atomic  oxygen 
number  density  only. 
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The  discharge  is  assumed  to  occur  in  a long  cylindrical  tube  when 
an  external  electric  field  is  applied  along  the  axis  of  the  tube.  The 
cylindrical  coordinates  are  (r,  ij),  2)  with  2 along  the  tube  axis. 

The  electric  field  in  the  tube  will  be = described  as 

E = [Er(r),  0,  Eq]  (2.4) 

The  average  velocities  of  the  three  species  are  3^,  q_,  and  The 
subscripts  represent  e',  0“,  and  O^,  respectively.  These  average 
velocities  will  be  modeled  with  the  following  components: 

Sg  = [Vg(r),o,UgJ ; q.  = [v_(r) ,o,u_] ; 3^=  [v^Cr),o,u^]  (2.5) 

The  axial  components  of  the  average  velocities  for  each  species 
are  u^,  u_,  and  u^.  They  are  the  drift  velocities  down  the  tube  and 
will  be  assumed  to  be  constant  or  at  least  that  their  derivatives  with 
respect  to  the  radial  distance,  r,  will  be  negligibly  small.  This  is 
a reasonable  assumption  for  a constant  potential  applied  across  the 
ends  of  a long,  narrow  discharge  tube. 

With  these  preliminaries  completed,  it  is  now  possible  to  examine 
the  three  species,  three  reaction  discharge  in  detail.  First,  the 
"general  discharge  equations"  will  be  reduced  to  "specific  discharge 
equations."  Second,  these  "specific  discharge  equations"  will  be 
partially  solved  to  obtain  the  "essential  discharge  equations."  Third, 
a magic  substitution  will  be  made  to  change  these  essential  differential 
equations  into  a tractable  set  of  algebraic  equations.  And  finally, 
this  set  of  algebraic  equations  will  be  solved.  Their  solution  will 
constitute  an  analytic  solution  to  the  three  species,  three  reaction 
oxygen  discharge. 
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Reduction  of  the  General  Discharge  Equations 


In  view  of  the  specific  setup  outlined  above,  the  general  dis- 
charge equations  consist  of  the  particle  conservation  equations,  which 
are: 


n V 

''e'^e  ^ "e''e  ~T~  " ' “^"e 


V n'  + n v'  + = oin„  - yn 

— • — — r p ' — 


v+n;  + n^v;  + — = 3ng 


(2.6) 

(2.7) 

(2.8) 


The  prime  indicates  differentiation  with  respect  to  the  radial 
distance  r.  The  radial  momentum  conservation  equations  are 

-1 


V v'  = (T.n  ) - - V V 

e e n m^'-  e e m r e e 
e e e 


vV 


r n ) ' - - V V 

- - ' m r - - 


'+''1  = n4:<V+^'  -"^r 


(2.9) 

(2.10) 
(2.11) 


The  axial  momentum  conservation  equations  are 


e e 

“ol 

(2.12) 

v_u^ 

m_  0 

- v_(u_  - 

“o) 

(2.13) 

v^u; 

= ^ 
m_^  0 

■ - 

%) 

(2.14) 

energy  conservation  equations  are 

Ve=-f' 

T CV  + 
e'-  e 

+ In  V 

r 3 e 

,(V^  * <^1) 

• ^ 

3 eo  e 

(2.15) 

I V 


J.  = - f T_(v.  + -ji^)  + ^.v.(v2  + u?)  - |;_^(T_-T^)  (2.16) 


9 


C2.17) 


And  Poisson's  equation  is 

= - He  - n_)  (2.18) 

0 

A detailed  derivation  of  results  which  lead  directly  to  these  equations 
is  contained  in  Appendix  A and  Appendix  B. 

The  primes  indicate  differentiation  with  respect  to  the  radial 
distance,  r.  The  quantities  n^,  n_,  and  n^  are  the  number  densities 
of  the  indicated  species.  T^,  the  Maxwellian  electron  temperature,  is 
in  units  of  ergs  and  therefore  is  normalized  by  the  Boltzmann  constant. 
T_  and  T^  are  the  energy  temperatures  of  the  indicated  species,  m^, 
m_,  and  m^  represent  the  masses  of  the  species,  q^  is  the  unit  elec- 
trostatic charge  and  is  the  permittivity  constant,  v_,  and 
are  the  collision  frequencies  of  the  indicated  species  with  the  neutral 
background  gas,  O2,  multiplied  by  the  fractional  momentum  loss  each 

species  suffers  per  collision.  Likewise,  v , v , and  are  the 

eo  -0  +0 

I 

collision  frequencies  of  the  species  with  OgS  multiplied  by  the  frac- 
tional energy  loss  each  species  suffers  on  impact. 

Now  let 


^e  ^ "e''e’  ^ V+ 

Where  j^,  j_,  and  are  the  radial  current  densities  for  each  species. 
The  particle  conservation  equations  become 

Jp/ 

j;  + = (B  - a)ng  + Yn_  (2.20) 

J:  + - an^  - yn.  (2.21) 


I 
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(2.22) 


j'  + ^+/ 

'‘+  ' I 


The  radial  momentum  equations  (2.9),  (2.10),  and  (2.11)  may  be 
simplified  by  assuming  that 


T n'  q 

V v'  « — ^ ® 
e e 


r . V V 

e e e 


V v'  « 


T^n'  q. 


’+  + 


n_m_’  m "r" 

T.n!  q 
-I 

m^“r' 


v.v:  « ::i=±, 


V V 


V V 
+ + 


(2.23) 

(2.24) 

(2.25) 


These  are  reasonable  assumptions  for  a discharge  tube  with  a 

1 cm  radius  operated  at  pressures  of  5 torr  or  less.  The  electron 

inequalities  may  be  demonstrated  using  representative  values  for  the 

parameters  involved.  For  electrons  with  a 1 eV  temperature,  T^  is  on 
-12 

the  order  of  10  erg.  Typical  radial  electron  velocities  are  about 

5 5-1 

10  cm/sec.  v^  for  a 1 cm  radius  tube  may  be  approximated  as  10  sec 

since  v^(o)  must  be  zero  to  preclude  non-singular  fields.  This  is 
probably  a generous  estimate.  It  is  assumed  that  electrons  will  re- 
combine at  the  tube  wall.  Then  n^/n^  may  be  approximated  by  unity. 

may  be  obtained  from  mobilities  and  is  approximately  10  sec  . 

-12  -27 

q^  is  1.6  X 10  erg/volt.  The  mass  of  the  electron  is  about  10  g. 

It  is  assumed  that  the  radial  electric  field,  E^,  is  much  greater  than 
-4 

10  volt/ cm.  Then  v^v^  is  much  less  than  the  quantities  on  the  right 

of  (2.23)  and  may  be  neglected  when  compared  to  them.  A similar 
demonstration  may  be  made  for  the  ions. 

It  is  assumed  that  at  equilibrium  T^,  T_,  and  T^  are  constant. 
Utilizing  (2.19)  and  these  assumptions,  the  radial  momentum  equations 
become 

11 


iTigV^^e  iv^e  r 


(2.26) 


m V - r 


(2.27) 


(2.28) 


It  has  already  been  assumed  that  u^,  u\  and  are  negligibly 
small.  Assume  that  the  axial  component  of  velocity  for  the  neutrals, 
u^,  may  be  neglected  when  compared  to  the  axial  components  of  velocity 
for  the  species,  because  only  the  charged  particles  are  accelerated 
by  the  applied  electric  field.  This  is  reasonable  for  a discharge  in 
which  the  fractional  ionization  is  extremely  small  or  several  orders 
of  magnitude  less  than  unity.  This  same  assumption  could  be  made  if 
the  neutral  particles  were  pumped  through  the  discharge  tube  with  a 
velocity  equal  and  opposite  to  the  net  velocity  imparted  to  them  by 
the  charged  particles.  Then  the  axial  momentum  equations  may  be 


wri tten 


*^6  ” '"e^g'-o 


(2.29) 


m V 0 


(2.30) 


y = - ^ 

+ 0 


(2.31) 


The  energy  equations  may  be  simplified  by  comparing  terms  as  was 


done  with  the  radial  momentum  equations.  Assume  that  u^  » v^,  which  is 
assuming  that  the  drift  speed  of  the  electrons  along  the  tube  is  much 


12 


greater  than  the  diffusion  speed  of  the  electrons  along  the  radii. 


Also,  assume  that  whereas  v «v  andu^  «v.  . This  is 

reasonable  since  O'  and  O2  ions  are  about  the  same  size  as  the  back- 
ground O2;  therefore,  0 and  O2  will  have  a large  energy  coupling  with 
the  O^.  The  electrons  are  very  weakly  coupled  to  the  O2.  Remembering 
that  Tg,  T_,  and  T^  are  constant  at  equilibrium,  the  energy  conservation 
equations  may  be  written  as 


T = 

e e 

eo 


(2.32) 


T = T 


(2.33) 


(2.34) 


Lee  obtains  this  result  by  a non-dimensional ization  of  the  same  general 
energy  equations  (Ref  9:1-7). 

The  original  "general  moment  equations"  have  been  reduced  to 
"specific  moment  equations"  for  the  three  species,  three  reaction 
oxygen  discharge.  These  "specific  moment  equations"  are  the  particle 
conservation  equations. 


jg  + = (e  - a)ng  + yn. 

j;  + = ctfig  - Yn_ 

j;  + = en^ 


(2.20) 


(2.21) 


(2.22) 


the  radial  momentum  conservation  equations 


7 if  * 


(2.26) 


= 


■''‘e  % 

-nl  - 


e e m v e r 

e e e e 


-T 


J = 


m V " m V -^r 


'^+  ~ ni_|_v^"+  iil^y^+^r 


the  axial  momentum  conservation  equations 


e m 0 
e e 


u = 


m V 0 


~ m^v^^o 


the  energy  equations 


*^e^e  2 

T..  = 

eve 

eo 


T-  = T , = 


(2.27) 

(2.28) 

(2.29) 

C2.30) 

(2.31) 

(2.32) 

(2.35) 


and  Poisson's  equation 

e‘  + ^ ^ (n^  - n^  - nj  (2.18) 

This  system  of  equations  is  used  by  Lee  at  the  outset  of  his  solution 
for  the  identical  discharge  (Ref  10:1). 

The  Essential  Discharge  Equations 

The  methodology  her®  will  closely  follow  Lee's  treatment  of  the 
three  moment,  three  reaction  gas  discharge  (Ref  10:1-8). 


, Ja.- 
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For  the  specific  discharge  system  above,  it  is  assumed  that  the 
background  gas  temperature,  T^,  is  known.  The  ionic  species  tempera- 
ture will  be  represented  as  T^. 

The  specific  discharge  system  of  equations  is  (2.20),  (2.21), 
(2.22),  (2.26),  (2.27),  (2.28),  (2.30),  (2.31),  (2.32),  and  (2.18). 
Thus  there  are  eleven  equations  in  the  eleven  unknowns  j^,  j , j^, 

Ug,  u_,  u^,  n^,  n_,  n^,  E^,  and  E^. 

Boundary  Conditions.  The  boundary  conditions  which  will  close 
the  system  must  be  specified.  In  order  to  retain  system  syimetry  and 
to  preclude  non-singular  fields,  it  is  necessary  that  at  r = 0 the 
following  be  true: 


ng(o)  = n;,(o)  = n;(o)  = 0 

(2.36) 

o 

II 

o 

+ 

II 

o 

1 

II 

o 

0) 

(2.37) 

E (o)  = 0 

1 

(2.38) 

Now  at  r = R,  the  wall  of  the  tube,  it  is  assumed  that 

of  the  charged  particles  occurs.  Also  the  net  current 

Thus, 

1 

recombination 

should  be  zero. 

, . j+(R)  - jg(R)  - j_(R)  = 0 

(2.39) 

n^(R)  = ng(R)  = n_(R)  = 0 

(2.40) 

This  system  of  eleven  specific  discharge  equations  (four  ordinary  linear 

» 

i differential  equations,  three  ordinary  non-linear  differential  equations, 

three  linear  algebraic  equations,  and  a single  quadratic  equation), 

^ eleven  unknowns,  and  eleven  boundary  conditions  constitute  the  model 
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for  the  three  species,  three  reaction  oxygen  discharge.  Only  seven 
boundary  conditions  are  required  for  a solution.  The  four  extra 
boundary  conditions  came  from  the  original  differential  equations 
(2.6)  through  (2.18). 


Beginning  the  Solution.  It  is  reasonable  to  assume  quasi - 
neutrality  in  the  positive  column  which  excludes  the  plasma  sheath 
near  the  wall  of  the  tube  (Ref  17:234).  Thus, 


n_^  = n_  + n^  (2.41) 

Only  the  following  equations  are  pertinent  for  beginning  a quasi- 
neutral solution: 


Jg  + '^e/^  = (6  - a)ng  + yn. 


y.  + -/p " 


-n 


~ m v^'^e  ” m v "e"r 
e e e e 


-T. 


J = 


0 -■ 


m V 


m V - r 


-T 


j.  = 


+'+ 


Now  let 


j = - 1 


Then  (2.20),  (2.21),  and  (2.22)  become 


(2.20) 

(2.21) 

(2.22) 

(2.26) 

(2.27) 

(2.28) 


(2.42) 
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(2.43) 


j + = 0 

Using  the  boundary  condition  j(R)  = 0 and  solving  yields 

j+  = j.  + O'e  (2.44) 

This  is  the  ambipolarity  result.  Also, 

j+  = J:  + J'e  (2.45) 

This  eliminates  equation  (2.22)  from  the  system,  and  n^  in 
equation  (2.28)  may  be  eliminated  by  using  the  results  (2.44)  and 
(2.41): 

(2.46) 

must  be  eliminated.  Equations  (2.26)  and  (2.27)  can  be 

written 


Now, 


-T„n'  m„v„ 
n E = ^ 6 _ e e • 

e r q q ■^e 

0 ^0 


-T  n ' m \) 
^0  ^0 


n n E = n n E 
- e r e - r 


(2.47) 

(2.48) 


(2.49) 


Therefore, 


V-^’e  ^ Ve"Je  " Ve"-  * 


(2.50) 
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If  equations  (2.26)  and  (2.27)  are  substituted  into  equation 
(2.46)  and  the  result  manipulated  slightly,  the  following  equation 
results: 

Je(V+  Ve^  ^ ^o^"e  " 

Then  the  entire  model  is  reduced  to  the  four  "essential  discharge 
equations": 

jg  + = (3  - a)ng  + yn,  (2.20) 

n:  + ^-/^  = on^  - yn.  (2.21) 

je("’+^+  + Ve^  * j-(V+  ^ " ■(’’’e  ^ ''’o^^e  " 

Ve^-^e  " Ve^-  ^ (2.50) 

The  boundary  conditions  require  ng(o)  and  n_(o)  to  be  bounded 
and 

n^(o)  = n;(o)  = 0 (2.52) 

ng(R)  = n_(R)  = 0 (2.53) 

This  is  the  basic  characteristic  value  problem  for  T^.  It  must 
be  remembered  that  the  production  rates  a and  3 depend  on  T^. 

The  Magic  Substitution 

The  essential  discharge  equations  consist  of  three  linear  first 
order  differential  equations  and  one  non-linear  first  order  differential 
equation.  Even  though  the  system  is  a basic  characteristic  value 
problem  for  T^,  finding  an  analytic  solution  to  it  appears  to  be 


insidiously  intractable.  Using  equations  identical  to  (2.20),  (2.21), 
(2.26),  (2.27),  and  (2.28),  Sabadil  does  find  a numerical  solution  to 
the  discharge  (Ref  17:236). 

However,  Lee  (Ref  10:5)  suggests  the  magic  substitutions 

^ (2.55) 

The  subscript  A indicates  all  species.  and  are  constant  current 
densities  and  number  densities,  respectively.  is  the  first  order 
Bessel  function  and  is  the  zeroth  order  Bessel  function.  The 
argument  of  these  Bessel  functions  is  defined  as 

C 

o I = (2.56) 

where  R is  the  radius  of  the  discharge  tube.  Also, 

Jo(c.)  = 0 (2.57) 

where  is  the  first  zero  of  J^.  This  is  similar  to  the  Schottky 
method  for  analyzing  a much  simpler  positive  column  (Ref  21:238-243). 

It  can  be  shown  that 

j;;  + ^X/r  = (2.58) 

nj^  = -a,N^J,(a,r)  (2.59) 

The  boundary  conditions  at  r = 0 and  r = R are  still  met. 

Substituting  (2.54),  (2.55),  (2.58),  and  (2.59)  into  the  essential 
discharge  equations,  and  assuming  that  T^  » T^  and  m^v^  » m^v^ 
results  in  the  system 
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(e  - a)Ng  + 


yn 


o.Jg  + 0 


= 0 (2.60) 


aN, 


yN  + 


a,J_  = 0 (2.61) 


^ “ V+'^e  ■ ^'"+^+  * ^ ° (2.62) 


-a.TNN  +a,TNN  +mvNJ  - 
'e-e  'oe-  ee-e 


m_v_NeJ_  = 0 (2.63) 


Thus,  Lee's  magic  substitution  has  changed  the  intractable  set  of 
linear  and  non-linear  differential  equations  into  a tractable  set  of 
three  linear  and  one  non-linear  algebraic  equations. 


Solution 

The  solution  to  the  three  species,  three  reaction  model  of  the 
oxygen  discharge  has  been  reduced  to  solving  a system  of  four  algebraic 
equations.  The  system  may  be  written 


Y 


-a  I 


0 


-v+ 


0 


-a 


-(m^v^  + 


(2.64) 


N 


N, 

Je 

N 

J 

‘'a 

e e 

- m v 
" 

< 

^e 

(2.65) 


N.  Op  '^- 

The  linear  system  (2.64)  can  be  solved  for  j;j-,  , and  Thus, 


3 

-a,Tg  + a,am_v_  + 
2a,  T^  + o(a,m_v_ 


(2.66) 


(2.67) 


t 


J 

N 


e 

e 


-2o,l^{a  + 6)  + SY(m+v+  + m v ) - a?T  y 


2a7T^  + yOiin  V 


J_  2a?TQa  - + a?TgY 


(2.68) 


A further  simplification  is  possible.  By  algebraically  carrying  out 
the  division  in  equation  (2.66),  it  is  seen  that 


+ - 2^' V] 

..  3_ 


(2.69) 


Then, 


N_ 

r = 

e 


(2.70) 


By  rearranging  equation  (2.67),  it  may  be  shown  that 


■’e.i 

O, 


(2.71) 


Results  (2.70)  and  (2.71)  may  be  substituted  into  equation  (2.65): 


e e 


•■j  ‘•2 


+ [a?(Te  - T^)  - mgVg(a  + B)  - ym.vj 


(2.72) 


This  is  the  characteristic  value  equation  for  Tg  . When  it  is  solved, 
the  three  species,  three  reaction  oxygen  discharge  model  that  has 
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been  postulated  will  be  determined.  Equation  (2.72)  is  solved  by 
iterating  it  with  T to  find  the  zeroes.  Once  the  equilibrium  T 
has  been  found,  it  is  possible  to  determine  the  remaining  unknowns 
in  the  model. 

Using  the  equilibrium  value  of  and  applying  equation  (2.32), 
the  axial  electron  drift  speed  is 


(2.73) 


Because  the  fractional  energy  loss  is  the  sum  of  the  elastic  as  well 
as  the  inelastic  processes,  determining  the  value  ^gQ  use  in 
this  equation  is  not  trivial. 

Equation  (2.29)  yields  the  equilibrium  axial  electric  field: 


^0  = 


eee 


Then  u_  and  u^  are  determined  by 


(2.74) 


u = 


m V 


(2.30) 


(2.31) 


Experimental  and  theoretical  information  on  the  equilibrium 
electric  field  versus  the  current  necessary  to  run  the  oxygen  discharge 
at  various  pR  values  is  available.  Dettmer,  in  his  thesis,  has  calcu- 
lated such  curves  for  the  oxygen  discharge  (Ref  2:38).  Of  course,  the 
model  of  the  discharge  herein,  with  the  resultant  equation  (2.74), 
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predicts  an  electric  field  independent  of  current.  Since  it  is 
possible  to  postulate  a value  for  the  discharge  current,  I,  at  the 
equilibrium  E^,  it  is  possible  to  compute  values  for  N^,  N , and 
(Ref  11): 


I = 27rq^  j (n^u^  - n_u_  - 
Jrs 


(2.75) 


Charge  quasi -neutrality  has  previously  been  assumed.  Utilizing  equations 
(2.54)  and  (2.55),  equation  (2.75)  becomes 

N 

I = 2TrqQNg[(u^  ‘ ^ F /rJ^(a,r)dr  (2.76) 

® Jn 


Then, 


2TTqoNeJ(u+  - Ug)  + ^u^  - u_)  R^J.(c,) 


(2.77) 


Equation  (2.77)  determines  and  thus  n^.  As  a result  N_  and  n_  are 

determined.  Then  n is  determined  by  charge  neutrality. 

Je 

The  ratios  and  ^ have  already  been  determined.  Given  Ng,  the 
e ' ,e 

radial  density  constants  Jg  and  J_  are  determined.  Therefore,  jg  and 
j_  are  determined,  and  results  from  applying  the  ambipolarity 
equation  (2.44). 

Finally,  the  radial  electric  field  E^  is  obtained  from  equation 
(2.47)  or  (2.48).  Using  the  magic  substitutions  as  well  as  (2.59)  it 


may  be  seen  that 


Ve 


(2.78) 


This  completes  the  determination  of  the  three  species,  three 
reaction  oxygen  discharge  model. 

The  Model  Inconsistency 

The  three  species,  three  reaction  discharge  model  which  has  been 
developed  is  actually  inconsistent.  But  this  is  the  normal  incon- 
sistency associated  with  the  Schottky  ambi polar  diffusion  model  with 
Bessel  profiles  (Ref  11).  Using  equation  (2.19)  and  the  magic  sub- 
stitutions (2.54)  and  (2.55),  it  is  easily  shown  that  the  radial 
velocities  of  each  species  are 


J. 

T 


o,r 


o,r 


(2.79) 


Equations  (2.78)  and  (2.79)  make  it  clear  that  and  v^  blow  up  at 
the  tube  wall,  r = R,  where  a,r  = c,,  the  first  zero  of  (see 
Figure  2-1 ) . 


Figure  2-1.  Inconsistencies  Due  to  Bessel  Profiles 
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Mitchner  and  Kruger  have  shown  that  the  ambipolar  diffusion  model 

begins  to  fail  at  a distance  from  the  tube  wall  which  usually  is 

negligibly  small  (Ref  13:150-154).  Lee  approximates  the  distance 
-4 

as  10  xR  (Ref  9:12),  The  region  near  the  tube  wall  is  the  plasma 
sheath;  and  modeling  it  into  the  discharge  would  require  a discharge 
model  much  more  sophisticated  than  the  simple  one  postulated  herein. 

It  is  assumed  that  the  distance  from  the  tube  wall  at  which  the  present 
model  begins  to  fail  is  very  small.  Therefore,  while  realizing  that  an 
inconsistency  does  in  fact  exist  in  this  model,  its  presence  will  be 
ignored. 

Summary 

It  has  been  demonstrated  that  an  analytic  solution  to  a simple 
three  species,  three  reaction  model  of  the  oxygen  discharge  does 
theoretically  exist.  The  agility  of  the  model  is  a result  of  Lee's 
magic  substitution  (Ref  10:5).  It  remains  to  be  demonstrated  that 
this- theoretical  solution  has  physical  significance. 

Earlier  it  was  pointed  out  that  the  dissociative  attachment  rate 
a and  the  total  ionization  rate  3 were  being  modeled  as  functions  of 
the  Maxwellian  electron  temperature  T^.  In  order  to  solve  the  charac- 
teristic value  equation  (2.72)  for  T^,  it  is  necessary  to  know  the 
functional  dependence  of  a and  6 on  T^. 
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III.  The  Dissociative  Attachment  Rate 


In  order  to  utilize  the  model  solution  obtained  in  Chapter  II, 

it  is  necessary  to  know  how  the  dissociative  attachment  rate  a depends 

on  the  Maxwellian  electron  temperature  T^.  The  dissociative  attachment 

reaction  makes  it  obvious  that  a is  highly  dependent  on  T : 

e 

O2  + e"  0"  + 0 (2.1) 

a as  a function  of  electron  energy,  E,  may  be  obtained  in  a 
manner  similar  to  the  method  employed  by  von  Engel  (Ref  21:292-293). 

A result  found  in  most  plasma  textbooks  and  demonstrated  by  Nielsen 
(Ref  15:26,  34)  is  that  the  probability  for  a dissociative  attachment 
event,  P^,  occurring  is 


Here  a(E)  is  the  total  collision  cross  section  for  dissociative  attach- 
ment as  a function  of  the  electron  kinetic  energy  E.  Assume  that  the 
electron  energy  distribution  is  the  Maxwellian 

N(E)dE  = ng  ^ dE  (3.2) 

Then  the  dissociative  attachment  rate  is 

o.(T,)  ■ /:  No^Pd(E)N(E)dE  (3.3) 

E^  is  the  threshold  energy  for  the  dissociative  attachment  process. 

N 

O2  is  the  background  O2  number  density.  Substituting  (3.1)  and 
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(3.2),  one  sees  that  the  rate  a satisfies 


First,  the  collision  cross  section  dependence  on  the  electron 
energy  will  be  determined.  Second,  using  the  determined  a(E),  the 
integral  of  equation  (3.4)  will  be  evaluated. 

Cross  Section  Dependence  on  Electron  Energy 

Cross  section  ^^(E)  versus  electron  energy  E data  for  dissociative 
attachment  in  oxygen  is  available  in  the  JILA  Report  (Ref  7:98-101). 
These  data  were  analyzed  with  a Hewlett-Packard  HP-67  calculator.  The 
data,  when  graphed,  assume  a bell-shaped  profile  (See  Figure  3-1). 

The  data  which  were  used  were  for  0^  at  a temperature  of  about  350°K. 
Data  for  0^  at  temperatures  as  high  as  1030“K  were  shifted  only  slightly 
toward  lower  energy  values. 

Following  von  Engel  (Ref  21:292-293),  ^^(E)  would  be  approximated 
as  the  best  linear  fit  to  the  leading  edge  of  the  data  curve.  In  most 
cases  this  is  quite  satisfactory,  because  it  is  only  the  meagre  supply 
of  electrons  in  the  tail  of  the  Maxwellian  electron  distribution  which 
make  it  past  the  threshold  energy  E^,  i.e.,  usually  the  tail  of  the 
electron  distribution  goes  to  zero  under  the  leading  edge  of  the  bell- 
shaped  aj(E)  curve. 

However,  in  this  model  it  was  not  clear  that  such  a linear  approxi- 
mation would  be  adequate;  therefore,  an  analytic  curve,  was 

determined  which  had  the  same  bell-shaped  profile  as  the  entire  data 
set.  Even  if  this  approach  were  needless  fine  tuning  for  a model  based 
on  a Maxwellian  distribution,  it  would  still  serve  to  limit  the  number 
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of  assumptions  made.  The  analytic  curve  was  determined  in  two  steps. 

First,  using  the  HP-67  "Curve  Fitting"  program  (Ref  6:03-01, 
L03-01)  and  adjusting  the  collision  data  to  the  program  coordinate 
system,  a reasonable  power  curve  fit  to  the  forward  edge  of  the  data 
was  found  to  be 

aj(E)  = 2.35  X 10'^^(E-3.8)^  (3.5) 

Where  3.8  eV  is  assumed  to  be  the  threshold  electron  energy  for  dis- 
sociative attachment,  E^. 

Second,  in  order  to  bend  the  curve  over,  o^(E)  was  assumed  to  be 
of  the  form 

aj(E)  = (3.6) 

2 -b 

Where  the  constants,  a and  c,  are  assumed  to  be  in  units  of  cm  ‘eV’ 

_2 

and  eV  , respectively.  This  equation  was  programmed  into  the  HP-67 
and  iterated,  varying  the  parameters  a,  b,  and  c.  The  best  curve  fit 
to  the  data  was  determined  by  eye.  It  was  found  to  be 

a^(E)  = (2.75  X lO"^^) (E-3.8)^e‘°-^^^^’^-®^^  (3.7) 

The  data  and  analytic  curve  (3.7)  are  shown  in  Figure  3-1. 
aCTe)/Mo 

Determination  of ^ 

The  dissociative  attachment  rate  “^^e^^N02  in  units  of 
3 -1 

cm  ‘sec  may  be  expressed,  with  the  aid  of  equations  (3.4)  and  (3.7), 
as 


28 


- f.', 


I 

i 

I 

I 

I 

I 


p 


1 


J 


.5 


♦ Schu;:!I209) 
GAjunai  (930) 

X \ X Papp(M6C) 

X A G Hende r 3on >384 1 1 


+ ♦ ♦♦♦♦♦♦ 
-j 1*  I 


ELECTRON  ENERGY  (aV) 


Figure  3-1.  Data  and  Analytic  Curve  for  Dissociative 
Attachment  Cross  Section  vs.  Electron 
Energy  in  Oxygen  (Ref  7:98) 


o.(Tg) 


2_  . -3/2 
fi  e 


/(2.75-10-'3)(E-3.8)V“-'S<E-3'8)^^^*\' 


In  order  to  express  all  energies  in  eV's,  it  was  necessary  to  convert 


E in  ergs  in  the  velocity  term  to  eV's: 


(E*)^(ergs)*^  = E^*(eV)*'"  x (1 .2657)(10'®) 


Combining  all  constants  and  using  (3.9),  equation  (3.8)  becomes 

so 

^ . (1-840411  X 1°'”)/E(E.3.8)  V'^"°e  -0-«(E-3.8)= 


Where  (1.840411  x 10~^^)  is  in  the  units  cm^*eV  ^‘sec’^  Now  let 


X = E-3.8 


(3.11) 


Then, 


(1. 840411x10"^^ ) f A -(0.19x^-h^) 

— e J ^ ^ ® 


No 


2 t/^2 


* (0.19x^-h^) 


2.  X 


/*0  vs/*»-//\»Yy 

'e  dx  (3.12) 


] 


The  evaluation  of  the  integrals  involved  in  equation  (3.12)  is  found 
in  Appendix  C.  Applying  the  results  of  Appendix  C to  equation  (3.12), 
one  may  show  that 


a{Tg) 

N02 


c, 


-3.8/T 

e 


erfc(^)  |c, 


^5  ^6 


+ c 


9 


•'ll  ^'12 

~2 V 

t:  Tf 

e e 


(3.13) 


An  HP-67  calculator  program  using  a series  sum  routine  was 
available  for  computing  the  erfc  function  (Ref  5:18-01,  L18-01).  This 
program  was  incorporated  into  a program  for  solving  equation  (3.13). 
However,  the  routine  overflowed  for  values  of  T^  £ 0.07  eV.  To  pre- 
clude this,  a small  T^  value  approximation  to  equation  (3.12)  was 
formulated  (Ref  11).  Equation  (3.12)  may  be  written 
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*3  8/ 

e ® /*  3 -0  19x^/ 

-jn;^  = c,^-^^7 / x'^e  x+3.8  e ®dx 

N02  ^ 3/2  ^ 


(3.14) 


For  small  values  of  the  e ^ term  will  dominate.  Thus, 


«(TJ  ? 3 -x/T 

-TITC^  = c,^--5-7 / 3.8x'^e  ®dx 

T 3/,  y 


■■2  T '^/2 
e 


(3.15) 


Replacing  e ^ with  1,  the  first  term  in  its  power  series  expansion, 
one  may  show 


a(Tg)  5/2  -3.8/Tg 

" '^13''’e  ® 


(3.16) 


The  HP-67  program  was  altered  to  include  equation  (3.16).  A 
“(T  ) , 

representative  list  of  ® 'N02  values  using  both  equation  (3.13)  and 
the  approximation  (3.16),  and  the  constants  c^  through  c.|3  are  in 
Sppendix  D. 

«(Te)/ 

It  is  of  interest  that  the  small  value  approximation  for  'N02 
is  about  one  order  of  magnitude  larger  than  its  exact  counterpart 
(3.13),  at  Tg  = 1.5  eV.  If,  in  the  final  solution  of  the  characteristic 
value  equation  (2.65),  it  is  found  that  the  equilibrium  T is  less  than 

c(rj  ® 

about  1.5  eV,  then  the  approximation  for  'N02  will  serve  to  ade- 
quately represent  the  dissociative  attachment  rate. 


Sunmary 


Using  the  JILA  Report  data  for  dissociative  attachment  in  oxygen. 


f 


an  analytic  representation  of  the  collision  cross  section  a^CE)  was 
modeled  as  equation  (3.7). 

This  o^(E)  was  used  to  determine  the  dissociative  attachment 

rate  which  is  modeled  by  equation  (3.13).  Equation  (3.16)  was  postulated 

to  provide  a low  T value  approximation  to  equation  (3.13). 

a(T  ), 

This  dissociative  attachment  rate  ® 'NO2  as  a function  of 
Tg  may  now  be  used  in  the  three  species,  three  reaction  oxygen  discharge 
equation  (2.72). 

However,  it  still  remains  to  determine  the  total  ionization  rate 

6 as  a function  of  T . 

e 


IV.  The  Total  Ionization  Rate 


Finding  the  functional  dependence  of  the  total  ionization  rate 
B on  the  Maxwellian  electron  temperature  T^,  is  the  final  analytic 
step  required  to  set  up  the  characteristic  value  problem  repre- 
sented by  equations  (2.64)  and  (2.72)  for  solution.  The  fact 
that  this  ionization  rate  is  also  highly  dependent  on  T^  is  indicated 
by  its  reaction 

O2  + e"  ^ O2  + 2e‘  (2.3) 

Determining  the  ionization  rate  B is  accomplished  in  exactly 
the  same  manner  as  determining  the  dissociative  attachment  rate  (Chapter 
III).  First,  an  analytic  representation  for  the  collision  cross  section 
0^. (E)  as  a function  of  the  electron  energy,  E,  is  determined  for 
ionization.  Then  this  result  is  used  in  equation  (3.4)  to  determine 
the  analytic  representation  for  the  ionization  rate. 

Cross  Section  Dependence  on  Electron  Energy 

The  data  used  to  determine  a^(E)  were  taken  from  the  JILA  Report 
(Ref  7:106).  The  threshold  energy  for  ionization  E^  is  12.063  eV. 
Because  of  this  relatively  high  E,  and  because  the  bell -shaped  data 
curve  does  not  begin  to  bend  over  until  the  electron  energy  is  over 
100  eV,  it  was  deemed  adequate  to  only  model  the  lower  portion  of  the 
forward  edge  of  the  data  curve.  The  data  points  are  in  Figure  4-1. 

For  comparison,  the  data  curve  for  dissociative  attachment  (Fig.  3-1) 
is  represented  on  Figure  4-1  by  the  small  triangle  to  the  left  of  the 
origin. 
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Figure  4-1.  Data  and  Analytic  Curve  for  Total  Ionization 
Cross  Sections  vs.  Electron  Energy  in  Oxygen 
(Ref  7:106) 

Using  the  HP-67  "Curve  Fitting"  program  (Ref  6:03-01,  L03-01)  to 
establish  an  initial  representation  of  a^. (E)  and  then  iterating  the 
constants  involved,  the  final  analytic  representation  was  found  to 
be 

a.(E)  = (4.0xl0'^^)e°-°^^^-  (6.8835145x10"^^)  (4.1) 

2 -1  17 

Where  o^-(E)  is  in  cm  , E is  in  eV"  , the  constants  (4.0x10"  ) and 

(6.8835145x10"^^)  in  cm^,  and  the  constant  0.045  in  eV“\  The 

analytic  curve  a^(E),  (4.1),  represents  the  data  well  up  to  about  30  eV. 

Figure  4-1  shows  the  data  and  the  analytic  curve  a^(E). 
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Determination  of 


Sf T ) 

The  ionization  rate  eVNo^  is  determined  in  exactly  the  same 
fashion  as  the  dissociative  attachment  rate.  Analogous  to  equation 
(3.8)  is 


■/[(‘'.0x10-'^)e‘’'“''=i(6.8835l45x10-”)l|j  1 
2 Ej  ■*  e 

(4.2) 

Collecting  the  constants  in  front  of  the  integral  and  rearranging,  one 
may  write  the  equation  as 


e(Tg) 


' T 


e 


Ee'^^e  ^^^®dE Ee  ^'^^^dE 

T 


(4.3) 


where 


1 

' a = 

2.6769x10'^ 

eV^^^  cm^  sec~^ 

(4.4) 

1 

1 

■ b = 

1 

1 

4.6067x10'^ 

eV"^^^  cm^  sec"^ 

(4.5) 

1 

i 

; c = 

i 

0.045 

eV"^ 

(4.6) 

* ■ = 

12.063 

eV 

(4.7) 

By  letting  x = E-E^.  and  carrying  out  the  integration,  it  may  be  shown 
that 


“Ej (y  "C) 

ae  ^ r 

NO2  j^3/2  L 


— E 

l-c]  (j-c)  ^ 

e JL  e J 


i/T„|V  4.C  T 


-bE 


(4.8) 


- ■ A- 


B(T  ), 

Some  representative  rates  for  ' NOo  are  tabulated  and  compared 
to  'N02  in  Appendix  D. 


Summary 

An  analytic  model  of  the  collision  cross  section  for  total 

ionization  as  a function  of  electron  energy  a^. (E)  has  been  postulated. 

This  collision  cross  section  has  been  used  to  determine  an  analytic 

representation  of  the  total  ionization  rate  as  a function  of  the 

S(T  )/ 

electron  temperature  T . Using  'No,  as  well  as  the  previously 

a(Te)/  ® ^ 

determined  NO2.  it  is  possible  to  compute  numerical  results  for 

the  three  species,  three  reaction  oxygen  discharge  that  was  modeled  as 
the  characteristic  value  problem  (2.64)  and  (2.72). 


i 

I 
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V.  Numerical  Results  for  the  Three  Species. 
Three  Reaction  Oxygen  Discharge 


The  characteristic  value  problem,  equations  (2.64)  and  (2.72), 

is  now  solved  by  using  the  results  from  Chapter  III  and  Chapter  IV 

which  provided  the  models  for  the  dissociative  attachment  rate  ^ ®^^No« 

6(TJ,  2 

and  for  the  ionization  rate  'N02. 

Method  of  Solution 

The  system  of  equations  used  to  compute  numerical  results  is 


i:: 

N 


2 2 
2a, T a - gyni.v  + a,T  y 

0 T + ^ 


e 2a, Tq  + YOim_v_ 


(2.68) 


ih  - « + ('] 

y T V 


(2.70) 


i = -L  . (:li) 

e e 


(2.71) 


°'Ve^N~^^  IP'(Te-To)-mgVg(a+6)-7m..'.3(^) 


- F(T,) 


(2.72) 


J N 

An  HP-67  program  was  written  to  evaluate  and  F(Tg). 

T was  adjusted  to  find  F(T  ) = 0.  The  program  was  initialized  by 

«(Tp)/ 

inputting  the  dissociative  attachment  rate,  ^ 'NOo,  and  the  total 

6<%)/  ^ 

ionization  rate,  'NOg. 
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F(Tg)  was  calculated  using  equation  (2.72)  for  a range  of  T^'s 
from  0.1  eV  to  20  eV  (see  Appendix  D for  a sample  list  of  values 
obtained  by  this  method).  By  observing  the  zero  crossings  made  by 
f^(Tg),  it  was  possible  to  establish  equilibrium  Tg  candidates.  Then 
P(Tg)  was  iterated  around  these  Tg  values  in  order  to  further  refine 
the  equilibrium  temperature. 

The  number  density  ratio  N_/Ng  was  used  as  the  criterion  for 
eliminating  non-physical  equilibrium  T values.  Since  the  ratio  N /N 

C — 0 

must  be  a positive  quantity,  any  equilibrium  Tg  values  which  had  an 
associated  negative  number  density  ratio  N /N^  were  discarded. 

In  order  numerically  to  solve  the  characteristic  value  problem, 
it  was  necessary  to  know  the  values  of  the  constants  a,,  T^,  nigVg, 
m_v_,  + . and  the  detachment  rate  y. 

Constants  Used  in  the  Solution 

The  final  solution  of  the  model  requires  that  the  values  of  the 
constants  in  equations  (2.68),  (2.70),  (2.71),  and  (2.72)  be  established. 

The  radius  of  the  discharge  tube  was  assumed  to  be  1.1  cm  (Ref  3). 
This  established  the  constant  in  the  Bessel  function  argument,  Jg(air) 
and  J, (o, r) , as 

o,  = 2.1862  cm"^  (5.1) 

The  background  gas  temperature,  T^,  was  assumed  to  be 

T = 300»K  (5.2) 

0 

In  order  to  establish  values  for  the  dissociative  attachment  rate 
a(Tg)  and  the  total  ionization  rate  6(Tg),  it  was  necessary  to  determine 
at  what  O2  background  gas  pressures  the  discharge  would  be  operated. 
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It  was  desired  tg  make  a comparison  with  Sabadil's  numerical  solution 
of  the  three  species,  three  reaction  oxygen  discharge;  therefore,  it 
was  decided  to  simulate  the  discharge  at  pressures  of  0.2  torr, 

0.5  torr,  1.0  torr  vRef  17:235-251).  5.0  torr  was  also  included. 

Sabadil  uses  electric  field  to  pressure  ratios.  E/p  in 
volt*cm"Vtorr,  of  40,  30,  and  20  for  the  operating  pressures  of 
0.2,  0.5,  and  1.0  torr,  respectively  (Ref  17:240).  Using  the  same 
E/p  values,  it  was  possible  to  determine  the  mobilities  for  the 
three  different  species,  with  the  mobility  relation 


McDaniel  and  Mason  give  a complete  set  of  mobility  vs.  E/p  curves 

for  the  ionic  species  in  background  0^  at  300°K  (Ref  12:292-293). 

These  curves  were  used  to  determine  the  values  of  m u and  m^v^  at 

all  the  operating  pressures  considered.  Sabadil's  electron  mobility 

values  for  electrons  in  0,  at  300'’K  were  used  to  determine  m for 

2 e e 

operating  pressures  of  0.2  torr,  0.5  torr,  and  1.0  torr  (Ref  17:240). 
Using  Brown's  data  for  the  electron  drift  speed  vs.  E/p  and  the  drift 
speed  mobility  relation  v^  = u^E,  it  was  possible  to  determine  the  value 
of  m^Vg  at  a pressure  of  5.0  torr  (Ref  1:57). 

The  final  constant  required  before  beginning  a numerical  calculation 
was  the  detachment  rate  y.  determined  by  the  reaction 

C - 0 - O2  + e’  (2.2) 

It  was  assumed  that  y was  equal  to  the  atomic  oxygen  number  density 
multiplied  by  an  appropriate  rate.  Garscadden  suggested  that  at 
equilibrium  the  atomic  oxygen,  0,  concentration  was  about  3%  of  the 
molecular  O2  concentration  (Ref  3).  An  appropriate  rate  for  the 
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reaction  was  taken  from  Dettmer's  dissertation  (Ref  2:263).  Thus,  | 

i 

Y = (0.03  No2)(1. 9x10'^ V-sec'^  (5.3)  I 

Table  VI  in  Appendix  D contains  a list  of  all  the  constants  used  J 

I 

for  the  solution  of  the  characteristic  value  problem  (2.72). 

Results  ' 

j 

The  results  from  the  solution  of  the  three  species,  three  reaction  j 

oxygen  discharge  model  are  displayed  in  Table  I. 

I 

Table  I 

Three  Species,  Three  Reaction  Model  Results 


Pressure 

in 

torr 

Equilibrium 

■'e 

eV 

N./Ne 

J./Ne 

in 

cm/ sec 

e e 
in 

cm/sec 

0.2 

3.67 

1.31  02 

-2.09  06 

3.14  06 

0.5 

2.21 

1.29  01 

-3.41  05 

5.22  05 

1.0 

1.50 

2.26 

-2.33  04 

4.25  04 

5.0 

0.99 

4.59  -01 

-5.18  02 

1.71  03 

J_/J 

-6.65 

-6.52 

-5.48 

-3.03 


e 

-01 

-01 

-01 

-01 


I The  reaction  rates  at  these  equilibrium  T^'s  are  displayed  in  Table  II. 

1 

The  negative  sign  on  the  ratio  J_/Ng  indicates  that  the  0~  ions  are 
! diffusing  radially  inward.  It  should  be  noted  that  if  one  assumes  that 

J^/Jq  is  approximately  equal  to  1,  the  radial  diffusion  speed  of  the 
electrons  is  given  by  Dividing  J /N„  by  N /N.,  the  approximate 

radial  diffusion  speed  for  the  O'  ions  may  be  found  (see  equation  (2.79). 


1 


Table  II 


Equilibrium  Reaction  Rates  for  the  Three 
Species,  Three  Reaction  Model 


Pressure 

in 

torr 

Equil ibrium 

^e 

eV 

ct(Tg) 

in 

sec'^ 

«Je) 

in 

-1 

sec 

Y 

in 

sec  ^ 

0.2 

3.67 

2.46  05 

2.31  06 

2.67  04 

0.5 

2.21 

4.03  05 

3.98  05 

9.17  04 

1.0 

1.50 

3.63  05 

4.20  04 

1 .83  05 

5.0 

0.99 

4.20  05 

2.61  03 

9.17  05 

Sutimary 

The  postulated  model  of  the  three  species,  three  reaction  oxygen 
discharge  has  produced  single  equilibrium  temperatures  which  are 
associated  with  physically  consistent  values  of  number  densities  and 
of  fluxes. 

Figure  5-1  indicates  that  the  electron  temperature  T^  at  equilibrium 
decreases  with  increasing  pressure.  It  may  perhaps  be  extended  to 
indicate  that  for  discharge  pressures  greater  than  5 torr  the  equi- 
librium electron  temperature  may  essentially  become  constant  around  1 eV. 

Figure  5-2  indicates  that  the  number  density  ratio  N_/Ng  also  de- 
creases with  increasing  pressure;  and  at  pressures  below  about  1 torr 
the  number  of  O'  ions  significantly  outnumbers  the  electrons. 

Figure  5-3  indicates  that  the  radial  diffusion  current  densities, 

J_  and  Jg  are  both  inversely  proportional  to  the  discharge  pressure, 
but  that  the  0”  current  density  decreases  faster  with  pressure  than  the 
electron  current  density. 
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In  Figure  5-4  the  curve  for  the  dissociative  attachment  rate  a, 
indicates  that  over  the  pressure  range  considered  the  production  of  ' 

0”  ions  is  saturated. 

The  0"ions  are  diffusing  radially  inward  toward  the  axis  of  the  ' 

discharge  tube.  This  result  concurs  with  that  of  Sabadil  (Ref  17:235). 

It  should  be  noted  that  the  method  of  solution  used  mobilities  to  , 

define  m_v_  and  These  mobilities  are  actually  functions  of  E/p. 

In  order  to  achieve  a closed  solution,  it  may  be  necessary  to  iterate 
the  entire  solution  based  on  the  equilibrium  electric  field  from 
equation  (2.72)  until  it  agrees  with  the  electric  field  used  to  deter- 
mine the  mobilities. 

The  primary  result  has  been  to  show  that  Lee's  method  for  attaining 
a solution  to  the  three  species,  three  reaction  oxygen  discharge  leads 
to  physically  consistent  results  in  some  cases  and  has  possibilities 
for  extension  to  more  than  three  species  (Ref  10:1-7). 
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in  cm  /sec 


Figure  5-3.  Equilibrium  Radial  Current  Density  Ratios  vs. 
Pressure  in  the  Three  Species  Model 
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Reaction  Rates  in  sec" 


Figure  5-4.  Equilibrium  Reaction  Rates  vs.  Pressure 

for  the  Three  Species,  Three  Reaction  Model 


VI . An  Equilibrium  Model  for  the  Four  Species. 
Four  Reaction  Oxygen  Discharge 


I 

1 

i 

i 


I 
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The  method  used  for  solving  the  three  species,  three  reaction 
oxygen  discharge  may  be  utilized  to  solve  the  four  species,  four 
reaction  oxygen  discharge.  The  new  species  is  the  negative  oxygen 
molecule  O2  . It  is  produced  by  the  charge  transfer  reaction. 

O2  + O’  ^ O'  + 0 (6.1) 

The  reaction  rate  associated  with  this  reaction  is  6. 

An  extra  reaction  and  ion  do  not  significantly  alter  the  method 
of  solution  previously  outlined.  Therefore,  only  the  cogent  differences 
will  be  outlined  here  before  displaying  the  form  of  the  characteristic 
value  problem. 


Solution 

The  specific  moment  equations  for  the  system  are  the  particle 
conservation  equations 


j / 

jg  = (e-ct)ng  + Yn_ 

J./ 

j'  + r = ouig  - (Y'''<5)n_ 
j;  + ^^/r  = 3ng 
32_  + ^^"^r  = 6n_ 

The  radial  momentum  conservation  equations  are 

■^e  ^ ifruT^e  ’ nuT^e^r 
e e e e 


(6.2) 

(6.3) 

(6.4) 

(6.5) 


(6.6) 
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-To  . 


m V - r 


(6.7) 


i = ^0  n"  + *^0  F 


(6.8) 


. 3 ~'*'o  ^ *^0  p 


(6.9) 


The  axial  momentum  equations  are 


u = — ^ 
e mgVg  0 


(6.10) 


m V 0 


(6.11) 


u = — 

+ "vT/o 


(6.12) 


2-  m2_V2_‘'o 


(6.13) 


The  energy  equation  is 


T - Ve^^e 
® '^eo 


(6.14) 


And  Poisson's  equation  is  used  to  close  the  system: 


0 


(6.15) 


..  ... 


It  is  assumed  that  the  0^  ions  are  at  the  same  temperature  as 
the  background  0^.  The  subscript  ^ identifies  parameters  associated 
with  the  0"  ion. 

To  eliminate  the  radial  electric  field,  E^,  the  following 
identities  may  be  used 


n n E = n n E 
- e r e - r 

(6.16) 

"2-"e^  ■ V2-^ 

(6.17) 

(6.18) 

The  system  can  be  reduced  to  the  six  equations 

J'e/ 

jg  + 'r  = (B-a)ng  + yr\_ 

(6.19) 

j'  + 'r  = an^  - (Y+6)n_ 

(6.20) 

^2-  * 

(6.21) 

jg(m^v^+mgVe)+j_(m^v^+m_v_)+j2_(m+u^+m2_V2.)  = ■’^e^Tg+T^j) 

■2V-  -^Va- 

" ^0^2-'’-  "’2-'^ 2-'^ 2-^-  (6.23) 

^o^'^-‘^"2-^"e'^e^eje^"-'^'^2-^  ^ TQ(n_+n2.) 

+ne(m_v_j_+iii2.V2-^2-^  (6.24) 

The  boundary  conditions  are 

49 


1/  ■ , 


(6.25) 


ng(0)  = n;(0)  = n'jO)  = 0 


ng(R)  = n_(R)  = n2_(R)  = 0 


(6.26) 


Using  an  extension  of  Lee's  magic  substitutions,  (2.54)  and 
(2.55),  the  system  becomes 


-0,  0 

0 

^2- 

N 

a - 0 

1 

'^e 

0 -Oi 

0 

l^e 

-a 

+N. 

y+  6 

0 0 

-a, 

J_ 

Ne 

0 

-6 

( m^v  ^-HiigV  g ) - ( m^v  ^+m_y 

_)-(m^v^+m2_V2_)| 

‘^2-1 

’ -'.(VTo) 

-2 '.To, 

plus  the  non-linear  equations 


(6.27) 


/^2-v  '^7. 

n7’^( it)"  '"2-'^ 2-^  n7^^n“ ^ 

e e e e 


(6.28) 


fN.  N2.'1  f J-  1 J J, 

[frg1i^p(V'^e^'^e^e(N;)J-'"-'^-<N7>-V^  = 0 (6-29) 


Equation  (6.27)  is  solved  for  the  ratios  tt-i  ti-t  '■"«  li — . 

\ \ Ng  * 

N_ 

in  terms  of  the  ratio  . Then  equation  (6.28)  is  solved  for  the  ratio 

e 

N. 

fj-.  Equation  (6.29)  is  the  characteristic  value  equation  for  computing 


N2.  Jg 


. jf-. 


the  equilibrium  T 's. 


2o?T„ 


Y 

a. 


-(:y+<s) 


a, 


= S(Tg) 


(6.32) 


It  has  been  assumed  that  T >>  T and  m v » m v . Then  the  solution 

e 0 - - e e 

of  the  four  species,  four  reaction  oxygen  discharge  is  the  following 
set  of  equations: 


(6.33) 


Where, 


Where, 


L “ ^ ( ^ R^) 


M = 'n_v.R^R3 


R,  = 


-T„  oon  V + 


2a?T, 


6-a 

a, 


(6.34) 

(6.35) 

(6.36) 

(6.37) 

(6.38) 
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1 

R,  = ^ 

3 a, 

1 

! 

(6.39) 

Snip  Vo  - m V (y+<5) 
s = -3  - 1 

2o?T„ 

(6.40) 

1 

S2  = ^ 

(6.41) 

1 s,  = -(ir+<s) 

1 3 a, 

I 

(6.42) 

. s = i- 

4 a, 

(6.43)  ' 

1 

Now, 

^■v4y 

± 

(6.44) 

n;  ■ "2  * s2[n^] 

1 

(6.45) 

1 

N;  = '^3  ^ ^31^] 

(6.46) 

1 

1 

(6.47) 

And  the  characteristic  value  equation  for  T is 

e 

r^-  ^2-ir  '^ei  j 

c c e ^ 0 0 

( 

(6.48) 

. 53 

t 

k.  ^ 

It  may  be  observed  from  equation  (6.33)  that  the  number  density 


ratio  TT"  two  values.  This  is  due  to  the  fact  that  r,-  is  the  root 

Ne  Ng 

of  the  quadratic  equation  (6,28).  Therefore,  the  characteristic  value 


equation  (6.48)  must  be  evaluated  for  both  the 


roots  for  a single  test  T . 

e 


N ■)  Tn  ■ 

and  the 

L'e  J - L e J + 


A program  was  written  for  the  HP-67  that  computed  F(T  ) for  both 

EN.1  fN.  1 No-  ^ h- 

jij-J  and  I J yielded  the  ratios  -jj— , and  Using 


this  program,  it  was  possible  to  begin  searching  for  the  equilibrium  T 


Numerical  Results 

In  order  to  solve  the  characteristic  value  problem  represented  by 
equation  (6.48),  it  was  necessary  to  determine  the  rate  6 for  the 
reaction 


O'  + O2  O2  + 0 


(6.1) 


Harrison  and  Moruzzi  have  investigated  this  reaction  rate,  6,  as  a 
function  of  E/p  (Ref  4:1244).  Their  results  were  used  for  interpolation 
to  the  lower  E/p  values  required  for  this  model,  i.e.,  40,  30,  20,  and 
10.  The  reaction  rate  5's  used  are  listed  in  Appendix  0 along  with  all 
of  the  other  constants  necessary  for  the  solution  of  equation  (6.48). 

The  results  of  the  computation  for  the  four  species,  four  reaction 
oxygen  discharge  model  are  listed  in  Table  III  and  Table  IV. 
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Table  III 

Four  Species,  Four  Reaction  Results 


Pressure 

in 

torr 

Equilibrium 

T in  e,, 
e V 

N_/Ng 

Nj/N, 

J /N 
^n® 
cm/sec 

J /N 
-k  e 
in 

cm/sec 

0.2 

5.51 

3.77 

5.48  03 

5.28  06 

7.72  01 

4.99  04 

0.5 

4.15 

5.25 

1.15  04 

4.38  06 

6.65  01 

6.91  04 

1.0 

3.54 

5.88 

2.12  04 

4.22  06 

3.28  01 

6.06  04 

5.0 

2.59 

5.11 

1.11  05 

4.52  06 

5.92 

8.11  04 

Table  IV 

Equilibrium  Reaction  Rates  for  the 
Four  Species,  Four  Reaction  Model 


Pressure 

in 

torr 

a 

^'’-1 

sec 

8 

in.i 

sec 

Y 

sec 

6 

in.i 

sec 

0.2 

5.51 

2.48  05 

1 .66  07 

3.67  04 

2.90  04 

0.5 

4.15 

6.33  05 

9.74  07 

9.17  04 

2.90  04 

1.0 

3.54 

1.21  06 

9.56  06 

1 .83  05 

2.25  04 

5.0 

2.59 

4.88  06 

1 .06  07 

9.17  05 

3.70  04 

I 

I 

1 

1 


• 

Again,  by  assuming  that  -r“  is  approximately  1,  the  radial  dif- 
fusion  speeds  of  the  electrons  and  of  the  negative  ions  may  be 
approximated  by  using  the  values  in  Table  III  and  equation  (2.79). 
The  results  are  displayed  graphically  on  the  following  pages. 


Summary 

Again,  the  postulated  model  has  produced  single  equilibrium  ' i 

temperatures,  T^.  They  were  regarded  as  physically  consistent  because 
neither  of  the  two  number  density  ratios,  N_/Ng  or  N2_/Ng,  was  negative. 

Secondarily,  these  T^  values  produce  other  physically  plausable  dis-  ^ 

charge  parameters.  I 

Figure  6-1  indicates  that  the  electron  temperature  decreases  with  i 

increasing  pressure.  This  may  be  due  to  the  fact  that  increasing  the  . 

pressure  reduces  the  radial  diffusion  to  the  tube  wall;  therefore,  the  ^ 

loss  rates  of  charged  species  is  reduced.  Then  the  available  electron 
energy  must  be  distributed  over  more  electrons,  and  the  average  electron 
temperature  thus  is  decreased.  Also,  if  one  assumes  that  the  electric  j 

field  varies  only  slowly  with  pressure,  the  electrons  have  shorter  paths  ; 

over  which  to  be  accelerated  at  the  higher  pressures. 

Figure  6-2  indicates  that  at  all  pressures  considered,  negative  i 

ions  outnumber  electrons.  But  now  it  appears  that  the  dominant  ion  is 
the  O2.  In  all  cases  it  far  exceeds  the  meagre  supply  of  O"  ions  in 
the  discharge.  The  number  of  0"  remains  fairly  constant  over  the 
pressures  considered,  but  the  number  of  0^  ions  increases  with  in- 
creasing pressure.  Thus,  it  may  be  concluded  that  the  charge  transfer 
reaction  (6.1),  which  produces  the  0^  ion,  is  the  dominate  reaction. 

This  conclusion  does  not  agree  with  the  experimental  observation  that 
the  0"  ion  is  more  numerous  than  the  0^  ion. 
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Figure  6-3  indicates  that  the  radial  diffusion  current  densities 
Jg/Ne^  J_/Ng,  and  relatively  constant  as  the  discharge 

pressure  is  increased.  Only  the  0"  ion  current  density  decreases 
slowly  with  increasing  pressure. 

Figure  6-4  indicates  that  as  the  pressure  increases  the  dissocia- 
tive attachment  rate,  a,  increases.  This  is  due  to  its  dependence  on  T. 

0 

and  the  neutral  species  density.  The  ionization  rate,  3,  initially 
decreases  but  then  remains  fairly  constant  over  the  range  of  pressures 
considered.  Since  the  detachment  rate,  y,  is  modeled  as  being  directly 
proportional  to  the  pressure,  it  increases  as  the  pressure  is  raised. 

The  charge  transfer  reaction  rate,  6,  remains  fairly  constant  as  the 
pressure  increases. 

In  this  model,  the  electrons  and  positive  and  negative  ions  are  all 
diffusing  outward  toward  the  tube  wall  under  the  conditions  postulated. 
This  indicates  that  the  outward  radial  flux  of  the  negative  particles 
balances  the  outward  radial  flux  of  the  positive  ions.  This  result  is 
supported  by  the  higher  production  rate  of  positive  ions  in  the  four 
species  model  than  in  the  three  species  model. 
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Figure  6-1.  Equilibrium  Electron  Temperature  vs.  Pressure  for 
the  Four  Species,  Four  Reaction  Model 
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Figure  6-2,  Equilibrium  Number  Density  Ratios  vs.  Pressure 
for  the  Four  Species,  Four  Reaction  Model 
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Figure  6-3.  Equilibrium  Radial  Current  Density  Ratios  vs. 

Pressure  for  the  Four  Species,  Four  Reaction 
Model 
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VII . Other  Models 


In  order  to  further  test  the  four  species,  four  reaction  model, 
some  other  reactions  were  added.  First,  the  three  body  attachment  re- 
action 

O2  + O2  + e~  O2  + O2 

was  included.  It  has  a reaction  rate  ic. 

Second,  the  charge  transfer  reaction 

O2  + 0 0"  + O2 

was  included  along  with  (7.1). 

Four  Species,  Five  Reaction  Model 

The  four  species,  five  reaction  oxygen  discharge  was  solved  in  an 
identical  manner  to  the  four  species,  four  reaction  model  discussed  in 
Chapter  VI.  The  equations  which  constitute  the  characteristic  value 
problem  are  found  in  Appendix  E and  the  constants  involved  in  its 
solution  are  in  Appendix  D. 

The  reaction  rate,  k,  was  taken  from  Dettmer's  dissertation 
(Ref  2:264).  It  was  represented  as 

K = (1.4xl0‘^^cm®*sec'^) [No2]^  sec"^  (7.3) 

In  this  case,  no  solutions  were  achieved.  It  was  possible  to  find  an 
equilibrium  electron  temperature,  e.g.,  at  1 torr  it  was  at  3.6  eV. 
However,  at  this  T^  the  number  density  ratio  N_/Ng  held  a steady  zero 
value.  This  was  considered  unacceptable.  Another  reaction  was  there- 
fore added. 
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(7-1) 


(7-2) 
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Four  Species,  Six  Reaction  Model 

In  an  attempt  to  eat  up  some  of  the  0^  ions  which  dominated  in 
the  four  and  five  reaction  models,  and  in  order  to  rejuvenate  the 
discharge  with  some  0”  ions,  the  charge  transfer  reaction  (7.2)  was 
added  to  the  system.  The  equations  which  define  this  characteristic 
value  problem  are  located  in  Appendix  E,  and  the  constants  involved 
are  in  Appendix  D. 

The  charge  transfer  reaction  rate,  e,  was  also  taken  from  Dettmer's 
dissertation  (Ref  2:263).  The  reaction  rate  was 

e=  (3.3x10  ^^cm^'sec"^  )(0.03)  [N02]  sec’^  (7.4) 

This  reaction  had  the  desired  effect  of  drastically  increasing 
the  population  of  0 ions,  but  it  also  disastrously  reduced  the  number 
of  O2  ions. 

The  radical  in  the  equation  analogous  to  equation  (6.33)  began 

! producing  complex  valued  ratios  for  the  number  density  ratio  N /N^  in 

1 the  regions  near  where  the  other  models'  equilibrium  T 's  were  located. 

e 

i The  present  treatment  of  the  four  species,  six  reaction  model  was 

i therefore  considered  unphysical  using  only  the  reactions  postulated. 

I 

I 
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VIII.  Discussion 


The  objective  of  this  thesis  was  to  obtain  reasonable  models  for 
the  three  and  four  species  oxygen  discharges  that  included  radial 
migration  of  the  electrons  and  ions.  Using  a three  moment  model  and 
making  appropriate  simplifying  assumptions,  it  has  been  shown  that 
both  the  three  and  four  species  oxygen  discharges  may  be  modeled  as 
basic  characteristic  value  problems  for  the  electron  temperature.  By 
using  a Bessel  type  substitution,  it  has  been  shown  that  the  differen- 
tial moment  equations  can  be  transformed  into  a simple  set  of  algebraic 
equations. 

The  facility  with  which  the  models  proposed  have  been  solved  is 
due  directly  to  the  analytic  technique  formulated  by  Lee  and  elaborated 
upon  in  this  thesis. 

It  has  not  been  claimed  that  the  solutions  obtained  are  unique. 

It  is  not  clear  at  this  time  that  there  are  no  other  solutions.  In 
fact,  a judicious  combination  of  Bessel  function  substitutions  is  the 
most  obvious  method  to  use  for  attempting  to  produce  other  solutions. 

The  models  analyzed  quantitatively  here  are  based  on  numerous 
assumptions  and  approximations.  This  is  undoubtedly  not  the  true  situ- 
ation, Dettmer  thinks  that  the  electron  distribution  is  between  a 
Maxwellian  and  a Druyvesteyn  (Ref  2:221).  Therefore,  it  is  not  sug- 
gested that  the  models  displayed  here  will  provide  quantitative  results 
for  the  actual  discharge. 

One  drawback  to  the  technique  used  to  solve  the  discharge  models 
is  that  it  cannot  handle  neutral  particles.  The  magic  substitution 
will  no  longer  smoothly  transition  the  system  of  differential  discharge 
equations  into  an  algebraic  set  of  equations.  However,  this  is  not 
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fatal  and  may  be  easily  handled  if  the  equilibrium  concentration  of 
the  neutrals  is  known. 

The  mathematical  approach  demonstrated  in  this  thesis  provides 
a rapid  and  a convenient  way  to  begin  the  analysis  of  gas  discharge 
models.  It  may  provide  useful  preliminary  results  in  parametric  studies 
which,  if  carried  out  with  more  accurate  models,  would  be  much  more 
expensive. 

The  three  species,  three  reaction  model  is  probably  better  than 
the  four  species  model,  because  it  corresponded  closer  to  what  has  been 
observed  in  the  actual  discharge.  Specifically,  the  0~  ion  and  not 
the  O2  ion  is  considered  to  be  the  dominant  species  (Ref  3). 

The  models  do  provide  an  analytic  tool  for  beginning  to  unravel 
the  trends  within  the  discharge.  Thus,  the  mathematical  technique  upon 
which  the  models  are  based  produces  solutions  which  indicate  that  the 
radial  profiles  are  Bessel.  Also,  the  models  indicate  that  ions  and 
not  electrons  are  the  numerically  dominant  particles  in  the  discharge 
at  low  pressure.  The  three  species  model  resulted  in  0"  ions  diffusing 
radially  toward  the  axis  of  the  tube.  Sabadil  arrived  at  this  same  re- 
sult in  his  numerical  solution  to  the  three  species  discharge  (Ref  17: 
235).  The  fact  that  the  equilibrium  electron  temperature  and  the 

N_ 

number  density  ratio,  decrease  with  increasing  pressure  is  a 

e 

physical  reality  (Ref  3).  Thus,  despite  its  simplicity  this  analytic 
technique  has  modeled  several  physical  observables. 

The  fact  that  the  four  species  model  began  to  apparently  diverge 
from  reality  is  probably  a result  of  using  too  few  of  the  thirty  odd 
reactions  In  which  0"  and  0^  participate  (Ref  2:263-266).  A more 
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appropriate  set  of  reactions  might  have  enlivened  the  four  species 
model.  Going  to  a five  or  six  species  system  is  perhaps  more  reasonable 
than  using  what  may  be  an  interim  four  species  system. 

To  this  end,  the  analytic  technique  which  has  been  used  to  model 
the  oxygen  discharge  is  perfectly  amenable  to  extension  to  systems  con- 
taining more  species  and  more  reactions  involving  them.  In  fact,  there 
is  theoretically  no  limit  to  the  number  of  species  that  could  be  main- 
tained in  a discharge  based  on  this  analytic  technique. 

Therefore,  the  outstanding  conclusion  to  be  drawn  is  that  Lee's 
magic  substitution,  (2.57)  and  (2.58),  and  straight-forward  general- 
izations of  that  substitution,  provide  an  extremely  simple  and  powerful 
tool  for  analyzing  multi -component  gas  discharges. 
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IX.  Recoirmendations 


Using  the  analytic  technique  developed  here,  the  four  species  model 
needs  to  be  improved.  It  may  be  useful  to  add  reactions  in  forward  and 
reverse  pairs. 

Throughout  the  entire  discussion,  the  Maxwellian  electron  tempera- 
ture dependence  of  m^v^  has  been  neglected.  The  models  were  perfectly 
amenable  to  entering  the  collision  frequency,  as 

Where  0^,1  is  the  average  momentum  transfer  cross  section  for  electrons 
in  oxygen.  In  order  to  do  this,  however,  it  is  necessary  to  know  the 
average  fractional  momentum  loss  per  collision.  In  this  same  vein,  the 
discharge  cannot  be  completely  solved  unless  the  energy  collision  fre- 
quency, is  known.  Neither  of  these  fractional  losses  were  known. 
Therefore,  it  is  recommended  that  a careful  analysis  be  done  to  determine 
Vg,  and  Vgg  as  functions  of  the  electron  temperature  Tg  in  the  oxygen 
discharge. 

Because  this  model  has  the  ability  to  theoretically  handle  any 
number  of  charged  particles  or  species,  an  analysis  which  includes  all 
the  species  in  the  oxygen  discharge  should  be  system! cal ly  attempted. 

Finally,  it  may  be  possible  to  combine  the  effects  of  gas  heating 
and  this  plasma  model.  Such  an  undertaking  would  represent  a synthesis 
of  the  thesis  done  by  Moore  and  the  models  contained  herein  (Ref  14). 
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Appendix  A 


Development  of  the  Moment  Equations  for  a Maxwellian  Gas 

Moment  equations  are  usually  the  cornerstone  in  the  development  of 
a plasma  discharge  model.  They  are  the  result  of  taking  successive 
moments  of  the  Boltzmann  transport  equation.  The  treatment  followed 
here  draws  upon  that  of  Rose  and  Clark  (Ref  16:112-126)  and  closely 
parallels  that  of  Lee  (Ref  8:1-7). 

Each  species  in  a plasma  discharge  has  its  own  Boltzmann  equation. 
The  derivation  of  the  Boltzmann  equation  is  found  in  most  plasma  physics 
textbooks.  Here,  equations  will  only  be  developed  for  one  species,  but 
similar  equations  may  be  written  for  other  species  in  the  same  dis- 
charge. The  Boltzmann  equation  may  be  written  in  non-relativistic  com- 
ponent form  as 


if  + V 

3t 


3f 
I 3x, 


3f  _ 


(A.l) 


Where  f is  the  distribution  function  for  the  species  involved  and  is  a 
function  of  x,  and  t the  position  of  a particle  in  position  space, 
the  velocity  of  a particle  in  velocity  space,  and  time,  respectively: 


f = f(x,  y_,  t)  (A. 2) 

The  parameters  x»,  v«,  v are  the  £th  components  of  position  and  the 
■c  4.  m 

£th  and  mth  components  of  velocity.  and  B^  are  the  Ith  and  nth 
components  of  the  electric  field  and  magnetic  field,  is  the 

permutation  symbol.  Z is  the  electrostatic  charge  number,  e.g.,  ±1, 
±2.  ^0  is  the  basic  electrostatic  charge.  The  mass  of  the  species 
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is  m.  The  term  on  the  right  hand  side  of  equation  (A.l)  is  the 
3f 

collision  term,  (— ) and  is  the  discontinuous  time  rate  of  change 

coll  ^ 

of  the  distribution  function  due  to  short  range  collisions. 

It  will  be  assumed  that  there  is  no  applied  field  and  that  the 
self-electric  field  is  negligible.  Equation  (A.l)  becomes 


+ Z 


9f 


^ ■ ^3t^ 


I 


coll 


(A.3) 


Particle  Conservation  Equation 

A particle  conservation  equation  results  when  the  zeroth  moment 
of  the  Boltzmann  equation  is  computed  by  integrating  (A.3)  over 
velocity  space: 


) dv  (A.4) 
coll 


The  number  density  for  the  species  is  defined  as 


n f d^ 


(A. 5) 


The  average  velocity  of  the  particle  is  defined  as 

a. 


= ly V f dv 


(A. 6) 


Now  it  is  assumed  that  the  distribution  function,  f,  is  piecewise 
continuous.  Then  (A.4)  may  be  written  as 


(|j)  dv 
V coll 


(A. 7) 
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Velocity  in  velocity  space  is  not  a function  of  position.  The  electric 
field,  E,  is  not  a function  of  velocity.  Thus,  using  (A. 5)  and  (A. 6), 
equation  (A. 7)  reduces  to  the  particle  conservation  equation: 


(A.8) 


Momentum  Transport  Equation 


The  momentum  transport  equation  is  obtained  by  multiplying  equation 
(A. 3)  by  mVj  and  integrating  over  velocity  space: 


= "1  f ' 

•V  ^ coll 


(A. 9) 


Which  becomes 


* srA 


'/'‘i 


(A. 10) 


Now,  using  equation  (A. 6)  and  the  fact  that  velocity  components  are 
independent  of  position  components,  (A. 10)  becomes 


I 


Define  a new  velocity,  U|^,  which  is  the  random  velocity  about  the  mean 
velocity,  qj^.  This  velocity  Uj^  is  called  the  "peculiar  velocity": 


% " "k  ■ \ 


The  second  term  on  the  left  of  (A.ll)  becomes 


(A. 12) 


Iv.v.fdv  + 

9xJ  j I - 3x, 


- |(Uj+q^.)(u^+q^)f  dv  (A. 13) 

jujU^fdv+^q^Ujfdy^+jqjU^fdv+qjq^jf  dv  (A. 14) 


_^f dv+nq^( q.  -q. ) +n( q^-q^) q ^ +fiq ^q^  ( A . 1 5) 


Therefore, 


^/vjv^fdv  - ^ 


Now,  the  kinetic  stress  tensor  appears  as 


(A. 16) 


UjU^fdv  =^Uj.u^f(x,  £+9^,  t)du  (A. 17) 


1 


Also,  the  third  term  on  the  left  side  of  (A.ll)  may  be  written 


I 3f  , _ ^'^Oc  3 

~m  m Z 3v^ 


IV  - dv 


(A. 18) 


(A. 19) 


Both  n and  are  independent  of  is  the  Kronecker  delta  function. 

Equation  (A. 19)  becomes 


!\  / V.  |f_dv  = - ^^E.n 


(A. 20) 


Now,  let  the  collision  term  on  the  right  side  of  equation  (A.ll)  be 


defined  as 


Jy  J coll 


dv  2 n A. 


(A.21) 


Then  using  (A. 16),  (A. 17),  and  (A.21),  equation  (A.ll)  may  be  written 


h ('"’ji  * 4 ["OjOi  * “j/]  * * " ''i 


»»5jnd1ng,  dividing  through  by  n,  and  using  equation  (A. 8),  equation 


OCCUMtS 


» • ’a" 


, 3o.,  Zq. 

t A, 

n ix-  m j J 


(A. 23) 
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It  is  assumed  that  the  distribution  function,  f,  is  the  Maxwellian 


f(x,v.t)  = n(2f^)3/2  e 


(A. 24) 


T is  the  Maxwellian  temperature  expressed  in  energy  units;  thus, 
embedded  in  T is  the  Boltzmann  constant.  The  quantities  n,  T,  and 
are  functions  of  and  The  kinetic  stress  tensor  may  be  simplified 
as  follows: 


=y  (v^.-q^)(Vj-qj.)f(x,v,t)d(v-g^) 

(A. 25) 

- — u^ 

"ij  ' ® ” ’‘S. 

(A. 26) 

a. . = <u.u .> 
ij  1 J 

(.A.27) 

The  quantity  <u.u.>  is  the  average  value  of  the  product  u.u.. 

Assume  that  u.  and  u.  are  uncorrelated.  Then  equation  (A. 27)  becomes 

* J 

O.J  = <u,>  <Uj>  (A.28) 

But  the  average  values  of  u.  or  u.  are  just  zero.  Therefore,  the 

^ J 

kinetic  stress  tensor  is  zero  unless  i is  equal  to  j.  This  result 
is  reasonable  in  light  of  the  fact  that  a Maxwellian  gas  will  not 
support  a shear;  thus,  only  the  diagonal  components  of  the  kinetic 
stress  tensor  are  non-zero.  These  three  components  may  be  evaluated. 


In  the  spherical  coordinate  system  (u,  4.,  0)  the  kinetic  stress  tensor 
becomes 


3 ^ y r r ■ 

n(^)  J J u^cose  sine  e du  de  d^)  (A. 29) 


Which,  when  carried  out,  is 


7 = ^ 

Kk  m 


(A. 30) 


Thus,  for  a Maxwellian  gas  the  kinetic  stress  tensor  is 


^ - nTf 

a . . = . . 

ij  m ij 


(A. 31) 


Then  the  momentum  equation  (A. 23)  becomes 


1 3 / nT \ 


29. 


3q. 


(A. 32) 


Energy  Conservation  Equation 

An  energy  equation  results  when  the  second  moment  of  the  Boltzmann 
equation  is  taken.  This  is  accomplished  by  multiplying  the  Boltzmann 
equation,  (A. 3),  by  and  then  integrating  over  velocity  space: 


f|''k''k''tlJ^l*  -5^ J' 


4 


v V dv 
-y^k''k3v^- 


coll 


(A. 33) 


f 
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Define  the  quantities 


ne  i f_^Vi,-q^)(V|,-qk)  f dv 


(A. 34) 

(A. 35) 


The  quantity  U is  the  average  kinetic  energy  per  particle.  Analogous 
to  the  peculiar  velocity  defined  by  equation  (A, 12),  the  quantity  0 
may  be  considered  a "peculiar  energy"  or  random  energy.  Expanding 
and  simplifying  (A. 35)  the  quantity  n0  may  be  written 


n0  = nU  - q^q^ 


(A. 36) 


Also,  utilizing  equations  (A. 17)  and  (A. 12)  in  CA.35),  it  is  seen  that 


/•v  _ m 

- I °kk 


(A. 37) 


Expanding  the  second  term  of  equation  (A. 33)  and  using  equations  (A. 12) 
and  (A. 17),  it  is  seen  that 


ml  3f  j..  _ m 3 

-j 


* "‘’k\‘i£ 


i 


"k“k“i  f dv  * q^o 


kk 


Now,  define  the  heat  flux  vector  as 


p _ m 
= I 


{wi 


f dv 


The  first  term  in  equation  (A. 33)  may  be  simplified: 


(A. 38) 


(A. 39) 


(A. 40) 
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Now,  define 

^ (A-41) 

Then  by  utilizing  equations  (A. 34),  (A. 38),  (A. 40),  and  (A. 41)  in 
equations  (A. 33),  it  is  found  that 

= Zq^E^nq^  * nP  (A.42) 

Now,  using  equations  (A. 36)  and  (A. 37),  equation  (A.42)  can  be 
rewri tten : 

= Zq^E^nq^  + nP  (A. 43) 

Expanding  this  equation  and  using  the  particle  conservation  equation 
(A. 8),  equation  (A. 43)  becomes 

IT  " ■ "^Icd  al^k  ■ ? *^k'’k^  2qQnq^E^+  nP  (A.44) 

The  result  of  contracting  the  momentum  transport  equation,  (A. 23),  with 


I 
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2 3V‘'j'’j>*2  sf 


Zq  E.q  . 

A.q. 

^ <Z  m j^j 


- ^ .q . 
J 


(A. 45) 


Using  the  result  of  (A. 45)  in  (A. 44)  one  finds  that 


3t  x9x^  n \(i  3x^^k  ” n 3x 


nr^)-n,Aj,j+P*  2 ^ (*-46) 


Equation  (A. 46)  is  the  energy  conservation  equation. 

Using  the  Maxwellian  distribution  results,  (A. 31),  equation 
(A. 37)  is 


0 = -J2-0  = JLr  ML)  = 1 

2n  kk  ^ m ^ Z 


{A.47) 


For  the  Maxwellian  distribution  the  peculiar  velocities  are  isotropic; 
thus,  equation  (A. 39)  for  the  heat  flux  vector  may  be  written 

(A.48) 

Equation  (A.48)  is  an  odd  function  of  u;  therefore. 


"I 


= 0 


(A. 49) 


Utilizing  {kAl)  and  (A. 49)  in  equation  (A. 46),  the  energy  conservation 
equation  is 


H 

3t 


+ 


2 

^mA.q . 
3 


Q m«  n 
n f’k^’k 


(A. 50) 
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Sunmary 

The  zeroth,  first,  and  second  moments  of  the  Boltzmann  transport 
equation  yield  respectively,  a particle  conservation  equation,  a 
momentum  conservation  equation,  and  an  energy  conservation  equation. 

If  it  is  assumed  that  there  is  no  applied  field,  that  the  self  electric 
field  is  negligible,  and  that  the  distribution  function  for  the  species 
is  Maxwellian,  then  the  following  three  moment  equations  result: 


' ^ 

1 3t 

\ 

( 


+ 


nq^)  = Q 


O 

Jt 


= _ 1 -i_ 

n 3Xj  ^ m ' 


Zq  E. 
^0  J 

m 


+ A. 


Qq. 


(A.8) 


(A. 21) 


ip 

Y 


Q m „ „ 
n 3 


S-T 

n 


(A. 50) 


For  a gas  containing  more  than  one  species,  these  equations  must 
be  written  for  each  species.  For  one  species  the  production  term  was 
represented  as  Q;  however,  for  a multi-species  system,  each  species 
would  have  a different  production  term.  Such  a system  of  equations  is 
related  through  their  production  terms. 


1 


80 


Appendix  B 

Equilibrium  Equations  for  a Two  Species  Gas  Discharqe 


Equilibrium  equations  for  a two  species  gas  discharge,  consisting 
of  electrons,  singly  charged  positive  ions,  and  a neutral  background 
gas,  are  developed  here.  The  system  will  be  assumed  to  have  cylindrical 
symmetry  as  shown  in  Figure  B-1 . 


Figure  B-1.  Cylindrical  System 
CRef  19:138) 


It  will  be  assumed  that  the  system  is  in  equilibrium.  It  will  also 
be  assumed  that  the  average  velocity,  has  only  radial  and  axial 
components,  which  are  functions  of  the  distance  r from  the  axis: 


a = [v(r) , 0,  u(r)] 


(B.l) 


It  is  assumed  that  the  electric  field  £ has  only  radial  and  axial 
components.  The  radial  component  is  a function  of  r,  but  the  axial 
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component  is  a constant: 


E = [E^(r),  0. 


(B.2) 


Particle  Conservation  Equations 

For  a two  component  system  in  equilibrium,  there  will  be  two 
particle  conservation  equations.  Using  equation  (A, 8),  the  electron 
conservation  equation  becomes 


Which  is 


^^(Oeq^)  = V.(n^)  = y = Q, 


n v'  + V n"  + = Q 

e e e e r e 


(B.3) 


(B.4) 


The  quantity,  Q^,  is  the  production  term  for  electrons.  The  quantity 
n^  is  the  number  density  for  electrons,  v^  is  the  radial  component  of 
the  average  velocity,  and  the  prime  indicates  differentiation  with 
respect  to  the  radial  distance  r.  A similar  equation  may  be  written 
for  ions: 


n.vT  + v.nr  + -3— = Q. 
11  11  r ^1 


(B.5) 


Momentum  Conservation  Equations 

In  equilibrium  the  momentum  conservation  equation  (A. 32)  is 

Tg  is  the  Maxwellian  electron  temperature,  and  m^  is  the  mass  of  the 
electron.  The  term  on  the  left  may  be  written  as  (Ref  11) 
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^ " ^(2  i'3.)  - g.  X V X 


(B.7) 


Now,  for  the  system  postulated 


1 

V(j  q-q)  = + e^Uj-jp 


(B.8) 


And 


Therefore, 


VX£=  f 


-gx(Vx£)  = 
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V 0 u 
e e 


/S  /N 

r 4, 


Ve  0 Ug 


0 - 


9r 


= -e 


4)  9r 


3u. 


(B.9) 


3u. 


= ®z''e^- Ve^ 


Then, 


= ®r''e''e  ^ ® ''e^^ 


(B.ll) 


The  first  term  on  the  right  of  equation  (B.6)  can  be  written  as  two 
component  equations.  The  radial  component  is  the  only  one  which  sur- 
vives in  the  discharge,  since  it  is  assumed  that  n and  T are  functions 

e e 

of  r only; 
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A similar  equation  may  be  written  for  the  ions: 

A^-^^  = - qj°^)  (B.15) 

The  superscripts,  (e),  (i),  and  (o)  represent  electrons,  ions,  and 
neutral  particles,  respectively,  is  the  collision  frequency  of 
electrons  with  neutrals  multiplied  by  the  fractional  momentum  loss 
that  an  electron  suffers  each  time  it  collides  with  a neutral  particle. 
The  fractional  momentum  loss  is  the  sum  of  the  elastic  and  inelastic 
momentum  losses.  Analogous  tov^,  the  termv^.  is  the  collision  fre- 
quency of  the  ions  with  neutrals  multiplied  by  the  fractional  momentum 
loss  suffered  by  an  ion  in  a collision  with  a neutral.  The  radial 
electron  and  ion  terms  are 


= -Vp(v^  - V ) (B.16) 


(B.17) 


V ) 
0 


It  will  be  assumed  that  the  neutral  background  gas  has  a 
negligible  radial  velocity  component  v^.  This  situation  could  be 
achieved  by  flowing  the  neutral  gas  through  the  discharge  tube  at  a 
proper  speed  dependent  on  the  radial  velocity  of  the  ions.  Thus, 
equations  (B.16)  and  (B.17)  become 


= -V  V 

r e e 

(B.18) 

r 11 

(B.19) 

The  axial  electron  and  ion  collision  terms  are 


(B.20) 

i(’>  = 

(B.21) 

The  last  term  in  equation  (B.6)  is  ‘ ^ j • This  term  will  be 
assumed  to  be  negligible.  This  assumption  means  that  newly  born 
particles  have  negligible  momentum. 

Using  the  results  (B.ll),  (B.12),  (B.13),  and  (B.16)  in  the 
equilibrium  momentum  conservation  equation  (B.6),  the  radial  momentum  con- 
servation equations  for  the  electrons  and  ions  become 


V V 

e e 


nVVe’ 

e e 


m. 


- V V 

e e 


(B.22) 


ViV.  + 


■<T,n,) 


n .m. ' 1 T 
1 1 


^''i 


(B.23) 
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Likewise,  the  axial  momentum  conservation  equations  for  the  electrons 


and  ions  become 


*^0^0 

V u = — ^ - V (u  -u  ) 
e e nig  e^  e o' 


ViU.  = 


2q 

0 0 


^-Vi(u.-Uo) 


(B.24) 


(B.25) 


Energy  Conservation  Equations 

At  equilibrium  the  energy  conservation  equation  for  electrons, 
(A. 50),  is 


3^e  2 P ? Qo  "’o  Q T 

^ • f e ^ ^ f ^ 


(B.26) 


Since  is  assumed  to  be  only  a function  of  r,  the  first  term  is 

3T 

(8-”) 


Also, 


It 

3 e 3x^ 


- Kv-a 


IUf  "e  ■ K < 'e  * ® ■ 28 ) 


Using  (B.18)  am  (B.20),  the  third  term  in  equation  (B.26)  becomes 
- K^j^j  ' - = |•e''e<''e  * “e  " %“o>  <8-2S) 


Since  u was  assumed  to  be  negligible,  certainly 


(v2  t u2)  » 


(B.30) 
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Therefore, 


! 


2 

3 


m A.q. 
e 


2 2 2 

Trfn  v^(v  + u ) 

3 0 0 ' g g ' 


(B.31) 


The  collision  term  P for  electrons  will  be  treated  in  a manner 
similar  to  Lee's  treatment  (Ref  7:8)  and  modeled  as 


P 


(e) 


-V 

eo 


(B.32) 


The  quantity  VgQ  is  the  collision  frequency  of  electrons  with 
neutrals  multiplied  by  the  fractional  energy  loss  suffered  by  the  elec- 
trons per  collision.  Also,  it  may  be  assumed  that  the  Maxwellian  elec- 
tron temperature  T^  is  much  greater  than  that  of  the  background  gas  T^. 
Therefore,  the  electron  collision  term  is 


P 


(e) 


= -V 


eo 


T 

e 


(B.33) 


Likewise,  the  collision  term  for  the  ions  may  be  written 

- V (B.34) 

Again,  is  the  collision  frequency  of  ions  colliding  with  neutrals 
multiplied  by  the  fractional  energy  loss  suffered  by  the  ions  per 
collision. 

The  last  two  terms  of  equation  (B.26)  may  be  written  as 

This  term  will  be  assumed  to  be  zero,  because  it  represents  the  approxi- 
mate difference  between  the  kinetic  and  thermal  energy  of  newly  bom 
particles.  This  difference  is  assumed  to  be  negligibly  small.  The  same 
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assumption  is  made  for  the  ions. 


Using  the  results  of  CB.27),  (B.28),  and  (B.31)  with  analogous 
equations  for  the  ions,  and  using  the  results  (B.33)  and  (B.34)  in 
the  energy  conservation  equation  (B.26),  the  following  two  energy 
conservation  equations  emerge  for  electrons  and  ions: 


V T'  = - (v'  + -r)  + u^)  - |v  T (B.36) 

ee  3e  e r'  3eee  e'  3eoe 


- iv,„(T,-T_,)  (B.37) 


Sutimary 

The  three  moment  equations  of  the  Boltzmann  transport  equation  were 
used  to  derive  a set  of  discharge  equations  for  a two  species  gas  dis- 
charge at  equilibrium.  The  moment  equations  assumed  only  that  there 
was  no  applied  ^ field,  that  the  self-electric  field  was  negligible,  and 
that  the  distribution  function  was  Maxwellian.  The  discharge  system 
had  cylindrical  symmetry,  but  none  of  the  vector  parameters  in  the  system 
had  an  azimuthal  component.  Also,  it  was  assumed  that  all  scalar  parame- 
ters were  functions  of  the  radial  distance  r only,  or  constants. 

The  particle  conservation  equations  are 


L 


n v'  + V n'  + 

= 

(B.4) 

e e e e r 

n Vi 

n.vr  + v.n:  + -V-J- 
11  11  i 

e 

(B.5) 

The  radial  momentum  conservation  equations 

are 
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'^e'^e  n m ^^e^e^ 
e e 

^0  r 

V V 

nig  e e 

(8.22) 

v.vT  = (T.n.) 

1 1 n.m^.  'll' 

+ ° ^ - V . V . 

m^.  1 1 

(B.23) 

The  axial  momentum  conservation  equations  are 


e e ni„ 


v.u.  = 


Zq  E 
0 0 


11  m.  1'  1 0 


The  energy  equations  are 


'e  X ^ 2. 


v. 


“o> 

(B.24) 

^o) 

(B.25) 

Ko% 

(B.36) 

(''X'  - l-io'^rTo) 

(B.37) 

Poisson’s  equation  may  be  used  to  close  the  above  system  of  eight 
differential  equations.  Thus, 

E'  + ^ = — (Z  n,-n^)  (B.38) 

r r r 1 e' 

0 

The  discharge  is  assumed  to  operate  at  pressures  of  5 torr  or 
less.  Therefore,  the  permittivity  of  free  space  is  almost  the  same 
as  the  discharge  permittivity. 
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i 

i 


It  should  be  noted  that  the  equations  derived  herein  are  not 

limited  to  two  charged  species.  They  may  be  extended  to  include  other 
ionic  species. 

I 

t>  , 

» 

i 


t 
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Appendix  C 


Evaluation  of  Dissociative  Attachment  Rate  Integrals 


In  Chapter  III  it  was  found  necessary  to  evaluate  the  integrals 


/ 

0 

a 

/ 


,4  ^-(0.19x2  * . p 


OB 

r. 

x"  e 


3 „-(0.19x2  * Vt^)  ^ 


(c.l) 


(C.2) 


These  integrals  may  be  evaluated  in  the  following  manner  (Ref  11). 
Let 


a = 0.19  ; b = y- 
e 


(C.3) 


By  completing  the  square,  it  may  be  shown  that 


dx 


(C.4) 


= ^-{/ax  + ^/2/a)2 


G = e / x-  e"^'“''  ' dx 


(C.5) 


Let 


w * /a  + ^^2/a  ; g * ^'^2/a 


(C.6) 


Then 


X . 

✓a 
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(C.7} 


Then  it  may  be  shown  that 


Now 


And 


“Y  (w^  - 4w^g  + 6w^g^  - 4wg^  + g^) 


{C.8) 


(C.9) 


(C.IO) 


•3  2 ^ 
3w  g + 


3wg 


(c.n) 


Let 


h~ 

9 


(C.12) 


And 


S “ 73/2  ^*"2  - 39^3  + 3g%  - g^Fg) 


(C.13) 


The  error  function  erf(k)  is  defined  as 


erf(k)  * — 

/if 


(C.14) 


And  the  complementary  error  function  erfc(k)  is  defined  as 


erfc{k)  = ds  = 1 - erf(k) 

^ k 


(C.15) 


Thus,  it  is  necessary  to  evaluate  the  five  integrals  F,  through  F-. 

I 5 

Using  Leibnitz's  rule  for  differentiating  an  integral  for  F^ , and 
it  may  be  shown  that 


^erfc(g) 


(C.16) 


=2  = *,  . ^ e-s'  * 1 e-a' 

■g 


(C.17) 


a> 

= J w^e’”^  dw  = ^ erfc(g)  + | e’^' 


(C.18) 


By  using  a differential  integrand,  it  may  be  shown  that 


o» 

= ywe-“ 


2 e-9‘ 

dw  = |- 


(C.19) 


And  it  is  easily  shown  that 


dw  . f 


erfc(g) 


(C.20) 


iifif 


Now  by  using  (C.16)  through  (C.20)  to  evaluate  and  it  may  be 
shown  that 

OB 


And 

G =y’xV(-'*'»'>dx . \ I [-  f - 9']  * ^ 


2 


(C.21) 


(C.22) 


When  the  parameter  g is  returned  to  its  original  form  using  (C.6)  and 
then  evaluated  using  (C.3),  it  may  be  shown  that 


OD 

pg-(0.19x2  + ^ ^ 


dx 


^ e erfc(^) 
e 


fd  + *^5  ^6  *^7 1 


. S . ‘‘lo  “'ll 
e e e 


Where  d^  through  d^^  are  constants. 


(C.23) 
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Appendix  D I 

Data  from  the  Three  and  Four  Species  Models  ' 


This  appendix  contains  numerical  data  used  in  the  solution  of  the 
following  models: 

3 Species  (e",  0",  Og);  3 Reaction  (a,  3,  y) 

4 Species  (e",  0~,  O^,  Og);  4 Reaction  (a,  3,  y,  6) 

4 Species  (e”,  0",  Og,  0^);  5 Reaction  (a,  3,  y.  <S,  <) 

4 Species  (e",  O",  Og,  Og);  6 Reaction  (a,  3,  y.  6.  <,  e) 

Parameters  Used  for  Evaluating  Equilibrium  T^'s 
in  the  Three  and  Four  Species  Models 
All  data  is  tabulated  in  Table  V. 


Table  V 

Data  Used  in  the  Models 


For  all  pressures: 


Tq  0 300-K  » 4.14132  -14  ergs 
R » 1.1  cm 
a,  “ 2.1862  cm“^ 


Symbol 

Units 

0.2  torr 

0.5  torr 

1 .0  torr 

5.0  torr 

E/p 

V/cm«torr 

40 

30 

20 

10 

NOg 

cm 

6.4371  15 

1.6093  16 

3.2185  16 

1.6093  17 

m v_ 

g/sec 

7.9307  -17 

2.1796  -16 

4.8693  -16 

2.9989  -15 

m^v^ 

g/sec 

1.8825  -16 

4.6031  -16 

8.8508  -16 

4.2842  -15 

1 

CM 

CM 

E 

g/sec 

1.7780  -16 

4.5641  -16 

9.4793  -16 

4.4428  -15 

m V 
e e 

g/sec 

3.6575  -18 

3.3168  -18 

2.9127  -18 

2.0025  -18 
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Table  V— continued 


I 

I 

I 

I 


Symbol 

Units 

0.2  torr 

0.5  torr 

1 .0  torr 

5.0  torr 

Y 

3.6691 

04 

9.1730 

04 

1 .8345 

05 

9.1730 

05 

6 

^1 

2.8967 

04 

2.8967 

04 

2.2529 

04 

4.0233 

05 

K 

5.8011 

02 

3.6258 

03 

1 .4502 

04 

3.6258 

05 

c 

mM 

6.3727 

04 

1.5932 

05 

3.1863 

05 

1.5932 

06 

Constants  through  ^*1® 

Dissociative  Attachment  Rate 

The  values  of  c^  to  c^^  ®'''®  Table  VI. 

Table  VI 
Constants  c-j  to 


Table  VI— continued 


Dissociative  Attachment  Rates  and  Ionization  Rates 


The  “^^e^/NOg  rates  computed  using  equation  (3.13)  and  the 
a(T  ) 

— ^approx  computed  using  equation  (3.16)  are  tabulated  in  Table  VII, 


6(TJ 


The 

Table  VII. 


N02  rates  computed  using  equation  (4.8)  are  tabulated  in 


Table  VII 

Dissociative  Attachment  and  Ionization  Rates 
as  a Function  of  Electron  Temperature 


Dissociative 

Attachment 

Rate  “^No«  in 

3 ^ 

cm  /sec 

Approximate 

Dissociative 

Attachment 

^NoT^  approx 

3 ^ 

in  cm  /sec 

Ionization 

Rate  ®^No,  in 

3 ^ 

cm'^/sec 

2.683649  -14 

3.712245  -14 

6.525468  -20 

2.687377  -12 

9.387083  -12 

1.767009  -14 

1.134319  -n 

9.180474  -11 

1.328730  -12 

2.131484  -11 

3.550301  -10 

1.258270  -11 

2.923942  -11 

9.069265  -10 

5.151058  -11 

3.452978  -11 

1.843090  -09 

1.379607  -10 

3.762679  -11 

3.247133  -09 

2.891278  -10 

3.913160  -11 

5.193019  -09 

5.189148  -10 

Table  VII--continued 


T 

e 

in 

eV 

Dissociative 

Attachment 

Rate  °^No-  in 

3 ^ 

cm  /sec 

Appro:' 'mate 
Dissociative 
Attachment 

^No^^  approx 

in  cm^/sec 

Ionization 

Rate  ^'^No,  in 

3 ^ 

cm  /sec 

IB 

3.954934  -11 

7.747175  -09 

8.389674  -10 

5.0 

3.925133  -11 

1.097011  -08 

1 .259933  -09 

5.5 

3.849647  -11 

1.491758  -08 

1.792519  -09 

6.0 

3.746108  -11 

1.964150  -08 

2.448331  -09 

6.5 

3.626376  -11 

2.519056  -08 

3.240551  -09 

7.0 

3.498369  -11 

3.161076  -08 

4.184566  -09 

Characteristic  Value  Computation 


Table  VIII  contains  a sample  listing  of  the  solutions  to  the 
characteristic  value  equation  (2.75)  for  a run  of  electron  temperatures. 
This  sample  is  for  the  1 torr  case  of  the  three  species,  three  reaction 

N./ 

model.  The  number  density  ratio,  'N^,  was  used  to  assess  whether  the 
equilibrium  T^  found  was  physical.  Notice  that  in  this  example,  the 
equilibrium  T^  is  found  between  1.4  and  1.6  eV. 

Table  VIII 

Characteristic  Value  Determination 


Tg  (eV) 

f (V 

N./Ne 

0.2 

-1 .260  -07 

-1.707  -02 

0.4 

-2.478  -07 

-3.343  -02 

0.6 

-3.735  -07 

-3.343  -02 

Table  VIII--continued 
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Characteristic  Value  Systems  for  the  Four 
Species.  Five  and  Six  Reaction  Models 


In  Chapter  VII  the  four  species,  five  and  six  reaction  models  were 
discussed.  Although  neither  system  was  solved,  the  setup  for  the 
characteristic  value  problem  is  included  here.  The  rate,  k,  is  good 
for  any  reaction  producing  an  O2  and  not  affecting  the  other  previous 
reactions.  The  rate,  e,  is  good  for  any  reaction  which  produces  an  O" 
ion  and  simultaneously  anihilates  an  0"  ion. 

Five  Reaction  Characteristic  Value  System 

Similar  to  the  four  species,  four  reaction  model,  the  five  reaction 
characteristic  value  system  is 

A,  ■ -‘-~Vl.^-4Kn' . Wl^-4KM  „„ 


’ 2K 

e 


Where, 

K - n'_v_S^S2  ~ "'2-2^4 
L = m_v_(R^S2  + 1135^)  - 
M = in_v.RiR3 


(E.2) 

(E.3) 

(E.4) 


Where 

_j  0(m_v_  + + <^2^2 

1 2T„  2afT„ 

0 0 

R,  . < 

2 a, 


(E.S) 

(E.6) 


1 

I 


I 
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,J'T  >». 


_ZJ£u. 


3 o I 


R = <_ 

'^4  0. 


- 1 


(E.7) 

(E.8) 

(E.9) 


S = 

2 o. 

S = ril.+<5) 

3 a, 


S 

^4  a. 


(E.IO) 

(E.ll) 

(E.12) 


Now, 


N N. 

^ = R]  + Si 

e 

J.  N_ 

“ ^2  ^2  ^N“^± 

e e 

J_  N_ 

= *^3  ^3 

e e 

J,_  N 

(n^)  “ *^4  ^4 


(E.14) 

(E.15) 

(E.16) 


And  the  characteristic  value  equation  for  is 
N N„  J Jo 


e e 


F(Tg)  » 0 


(E.17) 


Six  Reaction  Characteristic  Value  System 
As  before, 
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i 

i 


i 


■Vl2.4KM  ■^Vl^-AKM 

■ 2K  ’ ^fr>+  Ik 
e e 


(E.18) 


Where 


K = ni_v_S^S2  - '"2-^  2-^4 


L = m_v^(R^S2  + " '’’2-^2-*^4 


M = m V ^1 R3 


Where 


R.  = 


2 

ai(-o,Tg  + 3m^v^  + c<in_v_  + Km2_V2-) 


_ * 

2a7To  - <^1  + in_v_+  1^2 2-) 


R = 3 ~ g ~ < 


2 

o,(2Tga-Tge)  - e(2a-3)m^v^  + e(K-a)m2_V2. 

~ 3 

2a, - o,e(2m^v^  "’2-^2-^ 


«4  = 


of(2T^K-Tee)  + e(3-2K)m^v_j_  + e(a-K)m  v 
2 

2o,Tq  ^ a,e(2m^v^  + "'2-'^2-^ 


a,  [-2071^  - {Y+6)m_v_  + 6012 _V2_] 

S]  ~ 2 

2a, Tq  - a,e(2m^v^  + ni_v_  + m2_V2_) 


(E.19) 

(E.20) 

(E.21) 

(E.22) 

(E.23) 

(E.24) 

(E.25) 

(E.26) 


^3  = 


•2aiTQ(Y+6+e)  + 2e(Y+6)m^v^  +26an2_\^2_  YEn>2_V2. 

2a7Tjj  - o,e(2in^v^  + m_v_  + 012-'^ 2-^ 


(E.27) 

(E.28) 
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(E.29) 


2afTQ(6-e)  - 26e(m^v_^  + m_v_)  - Yan_v_ 

^4  ~ 5 

o.Tq  - a,e(2m^v_^  + m_v_  + m2_V2_) 


And 


J. 

W, 


uy 

+ 

II 

^nJ± 

e 

(E.30) 

N 

Vi 

e 

(E.31) 

N 

= "3  * ^3 

V- 

e 

(E.32) 

N 

= R,  + s, 

^ir^± 

e 

(E.33) 

And  the  characteristic  value  equation  for  T is 

e 


N.  J j j 

e ^ TT^  ^ ■"’2-'^  2-  ^ 

6 6 6 0 0 


F(TJ  = 0 (E.34) 
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